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Electroweak Sudakov logarithms at high energy, of the form (a/ sin 2 6w) n log" 1 s/M% w , are 
summed using effective theory (EFT) methods. The exponentiation of Sudakov logarithms and 
factorization is discussed in the EFT formalism. Radiative corrections are computed to scattering 
processes in the standard model involving an arbitrary number of external particles. The compu- 
tations include non-zero particle masses such as the t-quark mass, electroweak mixing effects which 
lead to unequal W and Z masses and a massless photon, and Higgs corrections proportional to 
the top quark Yukawa coupling. The structure of the radiative corrections, and which terms are 
summed by the EFT renormalization group is discussed in detail. The omitted terms are smaller 
than 1%. We give numerical results for the corrections to dijet production, dilepton production, ti 
production, and squark pair production. The purely electroweak corrections are significant — about 
15% at 1 TeV, increasing to 30% at 5 TeV, and they change both the scattering rate and angular 
distribution. The QCD corrections (which are well-known) are also computed with the EFT. They 
are much larger — about a factor of four at 1 TeV, increasing to a factor of thirty at 5 TeV. Mass 
effects are also significant; the qq — > tf rate is enchanced relative to the light-quark production rate 
by 40%. 
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I. INTRODUCTION 

Radiative corrections to high energy scattering pro- 
cesses have two powers of a large logarithm for each order 
in perturbation theory. These logarithms, referred to as 
Sudakov logarithms, lead to a breakdown of fixed order 
perturbation theory, and have to be summed to all orders. 
The Large Hadron Collider (LHC) has a center-of-mass 
energy of = 14 TeV, and will be able to measure 
collisions with a partonic center-of-mass energy of sev- 
eral TeV, more than an order of magnitude larger than 
the masses of the electroweak gauge bosons. Electroweak 
Sudakov corrections are not small at LHC energies, since 
a log 2 s/M^ z / '(4tt sin 2 9 W ) ~ 0.15 at y/s = 4 TeV. In 
this paper, we will apply effective theory methods de- 
veloped in two previous publications [l], to processes 
relevant for the LHC; in particular, we consider in detail 
dijet production, dilepton pair production, production, 
and squark pair production. In Rcfs. P, electroweak 
Sudakov corrections to the matrix element of an exter- 
nal current were found to be of order 10%. Electroweak 
corrections to LHC cross-sections are about four times 
larger. Naively, one factor of two arises because scatter- 
ing processes lead to four-particle operators, which have 



2 



(approximately) twice the radiative correction of the two- 
particle current operator. The other factor of two arises 
in squaring the amplitude to obtain the cross-section. 
Thus purely electroweak corrections at the LHC are sig- 
nificant, and resummed contributions must be properly 
included to obtain a reliable prediction for the cross- 
section. There are, of course, QCD corrections which 
are even larger, and are also included. 

There is an extensive literature on electroweak Su- 
dakov effects @, 3 i, 1, 0, i, 1 El, El El, El, EI El, El, 
E3, E^L The computations use infrared evolution equa- 
tions 5|, based on an analysis of the infrared structure 
of the perturbation theory amplitude and a factorization 
theorem for the Sudakov form factor [n| . These summa- 
tions have been checked against one-loop [H, EH, E3] ar± d 
two-loop [H, E3, El, El, E3 computations. 

The Sudakov logarithm log(s/M^, z ) can be thought of 
as an infrared logarithm in the electroweak theory, since 
it diverges as M\y,z —* 0. By using an effective field 
theory (EFT), these infrared logarithms in the original 
theory can be converted to ultraviolet logarithms in the 
effective theory, and summed using standard renormal- 
ization group techniques. The effective theory needed 
is soft-collinear effective theory (SCET) |H, HE [H HI , 
which has been used to study high energy processes in 
QCD [24[, and to perform Sudakov resummations aris- 
ing from radiative gluon corrections. 

This paper studies high energy electroweak corrections 
to processes relevant for the LHC, such as dijet pro- 
duction, dilepton pair production, tt production, and 
squark pair production, and expands on our previous 
works 0, 0, which will be referred to as CGKMl and 
CGKM2, respectively. In CGKMl we showed how to com- 
pute log s /Myy z corrections to the Sudakov form factor 
for massless fermions using EFT methods. In CGKM2 the 
results were generalized to massive fermions such as the 
top quark, including radiative corrections due to Higgs 
exchange. The corrections were computed without as- 
suming that the Higgs and electroweak gauge bosons 
were degenerate in mass. The Higgs corrections when 
expanded to fixed order agree with previous results of 
Mcllcs [18]. The electroweak corrections to processes in- 
volving four external particles are computed in this pa- 
per. We will show that the results can be obtained by 
summing the Sudakov form-factor results of CGKM2 over 
all pairs of external particles with appropriate group the- 
oretic factors. We also show how the results can be gen- 
eralized to processes involving an arbitrary number of 
external particles. 

There are different methods of counting the order of 
radiative corrections for the case of Sudakov corrections 
depending on whether one uses the amplitude or the log- 
arithm of the amplitude. We discuss this issue in detail 
in Sec. IIII1 where we also explain precisely which terms 
are summed in our computation. Roughly speaking, we 
use NLL running in QCD and LL running in the elec- 
troweak. The neglected terms are numerically less than 
1%. 



The paper is organized as follows: the outline of the 
calculation and notation is given in Sec. [H] The general 
structure of Sudakov double-logarithms, exponentiation, 
and the log-counting rules we use are given in Sec. IIIII 
We also discuss the numerical convergence of the pertur- 
bation series. The SCET formalism we use for our cal- 
culation is described in Sec. IIV1 including the formalism 
for Wilson lines needed in multi-particle processes com- 
puted using an analytic regulator [H, [26| . The calcula- 
tion of quark scattering and production is first calculated 
in a toy theory in section [V] Results are also given for 
massive quark production and squark production. The 
toy theory illustrates the theoretical tools needed for the 
standard model computation without the added compli- 
cations of a chiral gauge theory with three gauge groups 
and particles in many different gauge representations. It 
also illustrates how one can compute the radiative cor- 
rections for theories with scalar particles, such as super- 
symmetric extensions of the standard model. Some ob- 
servations on the factorization of amplitudes are made in 
Sec. IVII Radiative corrections in the standard model are 
given in Sec. IVIIl There are a total of eighty different 
amplitudes that are needed, which are computed in this 
section. Detailed numerical results and plots are given in 
Sec. IVIHl Appendix [A] discusses the box graphs needed 
for the high scale matching computation, as well as the 
crossing matrix needed for the case of identical particles. 
The parameter integrals we require in Sec. IVIIl are tab- 
ulated in Appendix [B] The top quark computation in 
CGKM2 was incorrect, and the corrected result is given in 
Appendix [Cl The numerical values change by about 1%. 



II. OUTLINE OF CALCULATION AND 
NOTATION 



The Sudakov logarithms are summed by integrating 
the renormalization group equations in SCET. The for- 
malism we use has been explained in detail in CGKM2. In 
this section, we outline the computation of four-particle 
processes; most of the results are well-known but will 
serve to define the notation we use in the rest of the pa- 
per. As in CGKM2, we first consider a toy gauge theory, a 
SU(2) spontaneously broken gauge theory with coupling 
constant a, where all gauge bosons have a common mass 
M. This is the theory used in many previous computa- 
tions [1, H, 0, H, H, E3] , and allows us to compare with 
previous results. The results will then be generalized to 
the realistic case of the standard model. When extend- 
ing the results of the toy theory to the standard model 
in Sec. IVIIl Higgs exchange effects will be included as in 
CGKM2. 

We consider two-to-two scattering at center-of-mass 
energies much larger than Mz- We will generically use 
Q » Mz to denote the energetic scale, and work in the 
regime where s,t, u are all of order Q 2 , so that one has 
hard scattering kinematics. Our results apply to high 
energy scattering processes at fixed angles, such as jet 
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(a) (b) 



FIG. 1: Tree level matching onto the intermediate effective 
theory. The full theory amplitude in (a) turns into scattering 
by a local operator in the effective theory, as shown in (b). 

production, but not to processes such as diffractive scat- 
tering. 

The scattering amplitude in the full theory arises from 
processes such as gauge boson exchange, as shown in 
Fig. The exchanged particle has virtuality of order 
Q 2 . At the scale /x ~ Q, we make a transition to SCET, 
which is an effective theory describing energetic particles 
with virtualities parameterically smaller than Q 2 . The 
full theory process is treated in SCET as scattering by a 
set of local operators, as shown in the right-hand graph 
in Fig. m 

i-Afuu = Q (/x) (p 4 p 2 | Oi (ji) \p!p 3 ) , (1) 

i 

where Oi are local SCET operators, and Ci (fi) are 
matching coefficients chosen so that the right-hand side 
reproduces the full theory amplitude up to power cor- 
rections of order M 2 /Q 2 . Power corrections can be sys- 
tematically included by keeping higher dimension oper- 
ators suppressed by powers of Q 2 . In our computation, 
we work to leading order in M 2 /Q 2 . The full and EFT 
have the same infrared physics but different ultraviolet 
behavior, and so we must introduce a set of matching 
coefficients, Ci(ji) which correct for the different short 
distance properties of the two theories. The matching 
coefficients Ci(/i) are computed by comparing on-shell 
matrix elements in the full and effective theories at a 
scale /i ~ Q. At this scale, infrared effects such as gauge 
boson and particle masses can be neglected, and so Ci (/i) 
can be computed using the unbroken gauge theory with 
massless particles. 

The coefficients Ci (/z) are evolved from /i ~ Q down to 
the scale /i ~ M using the SCET anomalous dimensions. 
The evolution equation for the matching coefficients in- 
volves the (matrix) anomalous dimension, 7^ , and is 

H^C i ^)= 7ij (fx)C j ( f x). (2) 

The anomalous dimension depends on the ultraviolet be- 
havior of SCET, and is independent of particle masses. 



Like the matching at Q, it can be computed using the 
unbroken theory with massless particles. In SCET, the 
anomalous dimension matrix can depend on log Q 2 /it 2 , 
so integrating Eq. |(3J) sums the Sudakov double loga- 
rithms. 

Once the coefficients Cj(/x) have been evolved down 
to a low scale of order M, we transition to a new effec- 
tive theory, which is also SCET, but with the massive 
gauge bosons integrated out. In our toy example, this 
new theory has no gauge interactions, since all the gauge 
bosons are massive. In the standard model, the transi- 
tion is from a theory with SU(S) x SU(2) x U(l) gauge 
bosons which we call SCETew to a new theory where 
the only gauge interactions are due to gluons and photons 
which we call SCET 7 . Operators Oi in SCETew are 
matched onto a set of operators Oi in SCET 7 . A single 
SU(3) x SU{2) x U(l) invariant operator Oi can break 
up into several operators Oi which are SU(3) x U(l) em 
invariant, but need not have full electroweak gauge in- 
variance. The SCETew ~~ > SCET 7 matching requires 
treating massive gauge bosons in SCET, using the for- 
malism developed in CGKM1, CGKM2. 

The operators in SCET 7 are evolved down to a scale 
set by the experimental observables of interest, and then 
used to compute the desired observables. For example, 
if one is interested in jet production, then the operators 
would be scaled down to fi of order the typical jet invari- 
ant mass. The operators can then be used to compute 
jet observables. This paper focuses on electroweak cor- 
rections, and we will not discuss this final step of the 
computation, since it is performed as discussed in earlier 
work [2t1 |. In our numerical results, we will choose this 
low energy scale to be 30 GeV. The electroweak correc- 
tions are not very sensitive to this scale, since the only 
effects below Mz are electromagnetic. The QCD cor- 
rections are scale dependent; the fi dependence in the 
SCET running cancels the /i dependence of the jet ma- 
trix elements to the order of the computation. We have 
not analyzed this in detail since we concentrate on elec- 
troweak effects in this paper. In Sec. IVIII1 only Fig. [T2l 
and Fig. 1131 have significant fi dependence. 

The bulk of the paper discusses the computation of the 
anomalous dimensions in SCETew an d SCET 7 , and the 
matching between SCETew and SCET 7 , which require 
SCET operators involving four-particles. We introduce 
the notation necessary to deal with an arbitrary number 
of particles. Most of the notation is standard to SCET, 
and we only discuss those features which are necessary 
for the extension to r-particles. 

The r energetic particles are described by SCET fields 
£ ni tPi labelled by momentum pi and light-cone direction 
rii, i = l,...,r. There are r light-cone directions rij, 
nf = 0, where — (l,n,), with a unit vector near 
the direction of motion of particle i. We will also define 
r light-cone directions rij by reversing the sign of space 
components of Uj, i.e. by applying parity to rij, = 
(1, — rij). Note that rij ■ m = 2. The momentum of any 
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FIG. 2: Vertex correction to the scattering amplitude in the 
full theory. 




gauge bosons denoted by A, with no labels, which couple 
to all the particles, analogous to the soft and ultrasoft 
fields introduced in NRQCD [28| . We work in the regime 
where the kinematic variables such as s,t are of order 
Q 2 , and the invariant masses of the final states are much 
smaller than Q 2 . The SCET power counting parameter 
is A = M/Q. The formalism is valid for observables that 
can be constructed out of variables in the effective the- 
ory, for which the reduction to effective theory vertices 
such as in Fig. © is valid. In particular, it is valid for 
jet observables and top decay observables at the LHC. 
Notation: We use the abbreviations 



-M 



log— 5-, L m = 
= l0g-5-, U 

M 



i 171 i 
log—, L Q 

M 

= log— , U 
A 1 



i Q 2 
:l0 V 

log^- 



U/t = log - = log(-s) - log(-i), 

U/u = log - = log(-t) - log(-u), 



i/i 



L ut/s2 = log — = log(-u) + log(-t) - 2 log(-s). (4) 



FIG. 3: Vertex correction in SCET. 



particle can be written as 

Pt = ^nf (fk • Pi) + ( n i ■ Pi) + Pi,±- ( 3 ) 

If THi is chosen to be exactly along the direction of pi, then 
pf L j_ = 0. The particles are energetic, with hi ■ pi ~ Q. 
In the case of only two energetic particles, one can work 
in the Breit frame where the particles are back-to-back, 
with hi = 7i2 and hi = m, so that one only deals with 
two null vectors n\ and hi, conventionally called n and 
h. 

Consider a radiative correction graph to the tree-level 
process Fig. [TJ such as the vertex correction shown in 
Fig. [2] in the full theory. The gauge boson exhanged 
between the two fermion lines still has virtuality of or- 
der Q 2 , and so the diagram behaves like the graph in 
Fig. [3J with the highly virtual gauge boson shrunk to a 
point. As is well-know, there are several different mo- 
mentum regions which contribute to the loop integral 
in Fig. O If the components of the gauge boson loop 
momentum are of order Q, then the gauge boson has 
virtuality of order Q 2 . This contribution is not present 
in SCET, and is included in the one-loop matching co- 
efficients at the scale Q. The other regions, which are 
included in SCET, are when the gauge boson is collincar 
to particle 1 (ni-collinear gauge boson), to particle 2 (n 2 - 
collinear gauge bosons), or is ultrasoft. The SCET the- 
ory thus contains nj-collinear gauge bosons for each par- 
ticle direction, i — 1, . . . , r, with momenta scaling like 
Pi, denoted by A nuPi with labels, as well as ultrasoft 



For scattering kinematics, s > 0, t < 0, and u < 0. 
All logarithms arise in the form log(— x — i0 + ) for x = 
s,t,u, so that log(— s — i0 + ) — logs — iir. Similarly, 
L s /i = log(-s) - log(-t) = log(-s/i) - m, and L t/s = 
log(— t) — log(— s) = log(— t/s) + in. This procedure can 
be used to find the branch cut of logarithms with negative 
argument which occur in the subsequent formulae. 



III. EXPONENTIATION AND LOG-COUNTING 

The exponentiation properties of Sudakov logarithms, 
and the relation between the renormalization group re- 
sults and those obtained by exponentiating fixed order 
computations was discussed in CGKM2. This section sum- 
marizes the results we need for our standard model cal- 
culation. 

The scattering amplitude A has an expansion of the 
form 1 

/ 1 \ 
cvL 2 cvL a 

a 2 L 4 a 2 L 3 a 2 L 2 a 2 L a 2 

a 3 L 6 



A = 



(5) 



1 For multi-particle scattering, A is actually a matrix of ampli- 
tudes, and matrix ordering is important. We discuss the simpler 
case of the Sudakov form factor, where A is a number. This is 
sufficient to study the exponentiating and log power-counting we 
need. The matrix case is discussed in Sec. IVII 
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where a represents a gauge coupling constant (a±, a-i 
or a s ), M is an electroweak gauge boson mass {M w or 
Mz) and Q ^> M is of order the center of mass energy 
of the scattering process, and L = logQ/M is the large 
logarithm. Each entry in Eq. (O has a numerical coeffi- 
cient, and the total amplitude is given by summing all the 
terms. The first row is the tree-level result, the second 
row is the one-loop contribution, etc. The a n contribtion 
has logarithms up to power L 2 ™. 

The logarithm of the scattering amplitude has an ex- 
pansion of the form [29|, H3, [M| 



log A 



/ aL 2 aL a 
a 2 L 3 a 2 L 2 a 2 L 



a 4 L 5 



(6) 



where the a n contribtion now has logarithms only up to 
power L n+1 , and the amplitude has been normalized so 
that its tree-level value is unity. The n th row can be 
computing using perturbation theory at n loops. There 
are far fewer coefficients in Eq. ([6]) than Eq. ([5]), so the 
form Eq. |6]) for log A is highly non-trivial. Equation (J6j> 
is referred to as the exponentiated form of the amplitude, 
since A is given by exponentiating the r.h.s. The first 
column gives the leading-log (LL) series, the second gives 
the next-to-leading-log (NLL) series, etc. 2 

The EFT computation naturally gives the scattering 
amplitude in exponentiated form. In general, there are 
several possible gauge invariants that contribute to the 
scattering amplitude, so that A is a matrix. The EFT 
computation gives the proper matrix ordering to be used 
for the exponentiated form of A. The difference between 
different matrix orderings can be computed using the 
Baker-Cambell-Hausdorff theorem. If X and Y are ma- 
trices, then 

e z = e x e Y 

Z = X + Y + [X,Y] + 1 [X, [X,Y\] + 1 % %X\) + 



where all the higher order terms are multiple commuta- 
tors of X and Y . If X and Y represent contribution to 
log A of the form Eq. ([6]), then X and Y are of order 
a n x,Y \_mx,Y wnere mx.Y < nx.Y + f • Thus one could in 
principle generate terms in log A of the form a n L m with 
m > n + 1 by reordering a matrix product using Eq. {7} . 
This does not occur, because, as discussed in Sec. EH the 
leading Sudakov series a n L n+1 is proportional to the unit 



matrix, and so drops out of the commutators in Eq. ||7J , 
so that the form Eq. © is preserved independent of the 
matrix ordering. 

When L is large, fixed order perturbation theory breaks 
down, and one needs to sum the logarithmically en- 
hanced higher order corrections. There are two inter- 
esting regimes relevant for the standard model, in which 
resummation is necessary. The first is the leading-log 
(LL) regime in which aL is of order unity. 3 This is the 
regime in TeV scale scattering for strong interaction cor- 
rections, where a — > a s . Using L ~ 1/a, the various 
terms in Eqs. ([5]) are of order 



/ 1 



A = 



a 

1 a a 2 



\ 



(8) 



/ 



Clearly the fixed order perturbation expansion breaks 
down, and higher order terms grow with inverse powers 
of a. To obtain a reliable value for the amplitude requires 
summing all terms along and below the diagonal, i.e. all 
terms of order unity or larger. The first super-diagonal 
gives the order a correction, the second super-diagonal 
gives the order a 2 correction, etc. 

The terms in the exponentiated form Eq. (J6j) are of 
order 



log A 



/—la 

a 

- 1 a a 2 

a 

- 1 a a 2 

a 
1_ 

a 

V i 



(9) 



The expression for log A has already achieved a partial 
summation of higher order terms. The largest terms are 
order 1/a, and there are no terms with higher powers of 
1/a. To obtain log A requires summing the first column 
(the LL series) and the second column (the NLL series). 
(7)The NNLL series gives order a corrections, the N 3 LL 
series gives the order a 2 corrections, and so on. While 
the NLL series is suppressed by one power of a relative to 
the LL series, it cannot be considered as a correction to 
the scattering amplitude A, since we have to exponentiate 
log A. If we write /„ for the N n LL contribution to log A, 
then 



log A = -/o + /i +af 2 + ... 
a 

= - [fo + afi + a 2 h + ■■■] 



(10) 



2 The LL, NLL, etc. counting used here is different from that used 
in fixed order calculations. The relation between the two is ex- 
plained in CGKM2, and after Eq. (1 1 1 ft in this section. 



3 Including loop factors of 4ir. 
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so that /i and ji are corrections to log A. However, 
1 



exp 



a 



-/o + fx + a/2 + 



= e° /o x e 71 x e Q/2 x 



(11) 



and exp /i can make a large change in A. Only /2 and 
higher can be considered as corrections to A. 

The counting discussed above is consistent with that 
used in renormalization group improved perturbation 



theory computations. In much of the literature, it is more 
common to use a different counting, which we denote by 
the subscript FO. The LLpo terms are those in A (not 
log A) of the form a n L 2n , the NLLpo terms are those 
in A of the form a n L 2n ~ 1 , and in general, the N fe LLpo 
terms are those in A of the form a n L 2n ~ k . In terms of 
fixed-order counting, Eq. |6j) can be written as 



log A = 



/ at 2 

a 2 L 3 

a 4 L 5 r 



" LLpo 

' NLLpo a 2 L 2 

N 2 LL FO a 3 L 3 
N 3 LL FO 



aL ~ NLLpo 

N 2 LL FO a 2 L 
N 3 LL F n a 3 L 2 



a ~ N 2 LL FO 
N 3 LL FO 



N 4 LL FO 



N 4 LL FO a 3 L~ N 5 LL F o a 3 ~ N 6 LL 



FO 



(12) 



and terms in A obtained by exponentiating are given by 
combining the powers of N. 

Note that with this counting, terms of a given series 
grow at higher order in perturbation theory, e.g. the 
N 3 LL FO terms are a 2 L, a 3 L 3 , a 4 L 5 , . . . , a n L 2n ~ 3 , which 
in the leading-log regime are of order a, 1, 1/a, 
l/a™~ 3 , and grow at higher orders. One can see this 
clearly from Eq. (fTTj) — fk+i is of order a k , and is small 
for k > 1, as are all terms in the expansion of expa fc /fc_|_i. 
However, the perturbation expansion for A contains the 
prefactor exp/ /a, and the terms (fo/a) n in the expan- 
sion of this prefactor for n > k mutiply the small terms 
in the expansion of e~xpa k fk+i to produce terms which 
are larger than unity, with a series of large contributions 
of alternating sign (since /o is negative). The problem 
is that the tree-level value A = 1 is not close to the true 
result for A; the leading contribution exp/o/a has an 
essential singularity at a = in the perturbation expan- 
sion. The second term exp/i also is not small. Only 
after these two contributions are factored out and prop- 
erly exponentiated does one have a reliable perturbation 
expansion. Summing all terms up to order N fc LLpo does 
not give a reliable calculation, because N fe+1 LLFo terms 
at order a r , r > k + 1 are larger than unity. It is essen- 
tial to properly exponentiate the /o and fx contributions 
to get a reliable expansion. Once this done, the higher 
order contributions are a small correction to the full am- 
plitude A. The amplitude A can be very different from 
the tree-level amplitude (a factor of 100 in our problem), 
and still be reliably computed in perturbation theory. 

The second regime we consider is the leading-log- 
squared (LL 2 ) regime in which aL 2 is of order unity. This 
is the regime in TeV scale scattering for electroweak cor- 
rections, with a — > ot\ t 2- Using L ~ 1/a 1 / 2 , the various 



terms in Eqs. ([5]) are of order 

( l \ 

1 a 1 ' 2 a 

A = 1 a 1 ! 2 a a 3 / 2 a 2 
1 

and in Eq. © are of order 



(13) 



log A = 



1 


a 1 / 2 


a 


a 1 ' 2 


a 


a 3 ' 2 


a 


a 3 / 2 


a 2 


a 2 







\ 



(14) 



/ 



The computation of A requires summing the first column, 
the Sudakov double- logs of order a™L 2 ". The remaining 
terms can be treated in a perturbative expansion. The 
second column gives the correction of order a 1 / 2 , the 
third column the order a correction, etc. The exponenti- 
ated form log A can be computed to order unity from the 
aL 2 term. The first correction, of order a 1 / 2 , is from the 
a 2 L and aL terms, the order a correction is from the 
a 3 L 4 , a 2 L 2 , and a terms, etc. We will refer to these as 
the LL 2 (leading- log-squared), NLL 2 , NNLL 2 , etc. con- 
tributions to log A. 

The scattering amplitude in the EFT computation has 
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the form 0, 
A = cxp 



2Ma(M))+X>i(a(M))log|L 



x exp ■ 



x expC(a(Q)) 



A(o(/z))log 



Q 2 



B(a(/i)) 



(15) 



Here exp C(a(Q)) is the high scale matching coefficient at 
Q 2 , 7 (/x) = A(a(^))log(/i 2 /Q 2 ) +B(a(/i)) is the SCET 
anomalous dimension between Q and M, exp D(a(M)), 
D(a(M)) = D (a(Af)) + L>i(a(Af))logQ 2 /A/ 2 is the low 
scale matching coefficient at M, a the gauge coupling 
constant (ai, ai or a s ), M is the electroweak gauge bo- 
son mass (Mw or M£) and Q » M is of order the cen- 
ter of mass energy of the scattering process. A is called 
the cusp anomalous dimension, and is linear in logQ to 
all orders in perturbation theory [H, • The low-scale 
matching expD has a single-log term D\ to all orders 
in perturbation theory fl], Q ■ The LL series is given by 
the one- loop cusp anomalous dimension, the NLL series 
by the two-loop cusp anomalous dimension, the one-loop 
value of B and the one-loop value of D\ , the NNLL series 
by the three- loop cusp, two-loop B and D\ , and one- loop 
.Do and C, and the N n LL series by the n + 1 loop cusp, 
the n-loop B and D%, and the n — 1 loop Dq and C. 
Eq. (fT5")) for the standard model, which we study in this 
paper, sums the QCD and electroweak corrections, in- 
cluding cross terms such as a s a\_2, o> s g 2 , or cx\^,9t which 
depend on mixed products of the Yukawa, strong and 
electroweak coupling constants. 



A. Absence of some terms in the Sudakov 
expansion 

In Eq. ([5]), we wrote the generic expansion for log A. 
In the standard model, one gets the form Eq. © where 
a n can be a product of the gauge or Yukawa couplings. 
It is interesting to note that not all possible terms are 
present. The leading Sudakov series in log A of the form 
a n L n+1 is given by integrating the one- loop cusp anoma- 
lous dimension with the leading order /3-function. The 
one-loop cusp anomalous dimension T(/i) is trivially a 
sum over the different gauge groups, since there can be 
no mixed terms like a s a\^ at one-loop, and because there 
is no Yukawa contribution to the cusp anomalous dimen- 
sion (see CGKMl). The one-loop gauge /3-function also 
does not mix different gauge couplings. Thus the leading 
Sudakov series is a sum of independent terms for each 
gauge group, with no mixed contributions, i.e. there are 
terms of the form a"L 2n , a"L 2 ™ and aV, L 2 ", but no terms 
of the form a"a^ 2 L 2n+2m for n,ra/0. 

The first contribution to the cusp anomalous dimen- 
sion which involves couplings from two-different gauge 
groups, and so cannot be written as the sum of contri- 
butions over individual groups, arises at four-loop or- 



der. 4 The two-loop (3- function also has contributions 
from two different gauge couplings. Thus at LL, the run- 
ning strong coupling a s only gets modified by terms of 
the form a s (a s L) n , but at NLL, one can have terms of 
the form a s (ai i 2a s L)(a s L)™(ai i2 L) m . The (ai^ot s L) fac- 
tor comes from one insertion of the two-loop /3-function 
in the renormalization group integration, and the other 
factors come from using the leading-order /3-functions for 
the remaining integration of a s and 0:1,2- 

Using the above, and noting that the matching con- 
ditions and non-cusp anomalous dimensions have all al- 
lowed terms, one finds that one can get all possible terms 
in Eq. ^ for the N 2 LL and higher series (third column 
and beyond). For the LL series (first column), all terms 
have a single gauge coupling. For the NLL (second col- 
umn), all terms can occur with the exception of the o 2 L 2 
contribution, which can only have a single gauge cou- 
pling, so that terms such as a s a\L 2 are absent. 



B. Terms included in the computation 

In the standard model, the radiative corrections in- 
volve the strong coupling a s and the electroweak cou- 
plings ai.,2 • For log-counting purposes, we assume that 
the strong coupling is in the leading-log regime, and 
the electroweak couplings are in the leading-log-squared 
regime. Let a be the log-counting parameter. Then 
L ~ 1/a, ai.2 ~ a 2 , a s ~ a. The top-quark Yukawa 
coupling is also treated as the same order as the elec- 
troweak couplings, g 2 ~ ai^ ~ a 2 - The terms in log A 
are given in Eq. (JSJ), but now each a can be either a 
strong coupling a s of order a or an electroweak coupling 
«i.2 of order a 2 . The order of terms with all couplings 
equal to a s are given by Eq. ^ with a — * a, those with 
one coupling 0*1,2 and the rest a s are given by Eq. © 
with a x (a — > a), etc. The leading terms of order 1/a 
in log A are given by summing the a™L n+1 terms, i.e. 
the leading-log QCD series. The order 1 terms are given 
by summing the a™L n and ai ! 2a™ _1 L™ +1 terms, i.e. the 
NLL QCD series and the LL series with one power of the 
electroweak coupling. The order a corrections are given 
by summing the a™L n_1 , ai i2 Os _1 L n and of 2 o™ _2 L™ +1 
terms, etc. In the exponentiated form Eq. ([6]), one only 
needs to include electroweak corrections at low orders, 
so that summing terms to order unity only require one- 
loop electroweak computations, to order a only requires 
two-loop electroweak corrections, etc. In contrast, the un- 
exponentiated form Eq. ([5]) of fixed-order computations 
requires electroweak corrections of arbitrarily high order 
to sum all terms of order unity or larger. 

In the numerical results of Sec. IVIII1 we include the 



4 i.e., r(o! s , 0:1,02) = r s (o s ) + ri(afi) + r 2 («2) up to three-loop 
order. We would like to thank Z. Bern and L. Dixon for helpful 
correspondence on this point. 
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one-loop QCD, electroweak and Higgs corrections, as well 
as the two-loop QCD anomalous dimension [34[ and two- 
loop running of the gauge coupling constants. This in- 
cludes the entire one-loop correction to the scattering 
amplitude, as well as all higher order corrections which 
are formally of order l/a or a? . The terms we neglect are 
order a or higher in the log-counting, and at least second 
order in the gauge couplings constants a s ,i,2- The error 
due to the neglected terms is numerically less than 1% in 
the rate. 

In terms of the commonly used fixed order counting, 
we have included all LLpo and NLLpo terms for both 
the QCD and electroweak corrections. In addition we 
have included all NNLLfo of the form a"L 2n ~ 2 and 
Using the counting that a s ~ a and 



I 2n-2 

Q1.2L 



«1.2 



, , and counting a n |_ - as N LLpo, we have 
summed all terms of order N 3 LLpo- hi terms of the ex- 
ponentiated form Eq. ©, which is the form given by 
SCET and used for the numerics, we have included 



/V V 



log A 



V 



\ 



\J not a\ 2 L 2 only 
y/ not a\ 2 L 3 x 
x/ not af 2 L 4 



V 



(16) 



/ 



where \J means all terms have been included, x means 
no terms have been included. The largest terms omit- 
ted are a\ 2 L 2 , a 3 L 2 and a s a\^, and are estimated to 
be (a/(7rsm 2 6W) 2 L 2 ~ 0.006, (a s /7r) 3 L 2 - 0.003 and 
a s a/{ir 2 sin 2 W )L ~ 0.003 using L - log(4TeV) 2 /M| ~ 
7. This gives a sub-1% error. The a\ 2 L 2 term arises from 
the two-loop electroweak cusp anomalous dimension, and 
the a 3 L 2 term from the three-loop QCD cusp anomalous 
dimension. These are known, and could be easily in- 
cluded in the computation. We have checked that these 
change the rates by less than 1%. 



IV. SCET FORMALISM AND WILSON LINES 

In SCET, m collinear gauge bosons can interact with 
particle 1, or with the other particles in the process. The 
coupling of rii-collinear gauge bosons to particle 1 is in- 
cluded explicitly in the SCET Lagrangian. The particle- 
gauge interactions are identical to those in the full-theory, 
and there is no simplification on making the transition to 
SCET. However, if an ni-collinear gauge boson interacts 
with a particle other than 1 (pick particle 2 for definitc- 
ness), then particle 2 becomes off-shell by an amount of 
order Q, and the intermediate particle 2 propagators can 
integrated out, giving a Wilson line interaction in SCET. 
The form of these operators was derived in Ref. [22j • We 
will use the definitions 



W(«2) 



exp 



(17) 



which is the expression given in Ref. [13] with the replace- 
ment n — ► m, n — * ri2' The gauge generators T A are in 
the representation 9^2 of particle 2. The subscript n\ 
is a reminder that the Wilson line contains ni-collinear 
gauge fields, and the superscript (n 2 ) is a reminder that 
the integration path is directed along 712, and that the 
gauge generators are in the representation of particle 2. 



("2 



is a e?2 x 0I2 matrix where di is the dimension 



of 9^2, and transforms under ni-collinear gauge transfor- 
mations as 



ur(n 2 ) 

* r ni 



ah 



ac 



(18) 



where is the gauge transformation matrix in the 9^2 
representation. One can similarly define Wn™ 3 ^ for any 
pair ij of particles, with i ^ j. It is convenient to treat 
all gauge indices as incoming, i.e. an outgoing fermion 
line in the gauge representation £H will be treated as an 
incoming fermion in the representation SH. 

A generic gauge invariant local operator in the full 
theory can be written as the gauge invariant product of 
fields, 



(19) 



where \i,a t is thi,a z for incoming particles, Xi,a t = ^], a< f° r 
outgoing particles, and c is a Clebsch-Gordan coefficient. 
Xi transforms as EHi for incoming particles, and as 91^ 
for outgoing particles. The indices on are gauge indices, 
and c {{a-i}) = c{a\, . . . , a r ) is the Clebsch-Gordan coef- 
ficient for combining the product of fields into a gauge 
singlet. For ni-collinear gauge couplings, the field xi m 
Eq. ([19]) can be replaced by the SCET field and 
the other fields are replaced by Wilson lines. Collinear 
gauge invariance implies that the operator Eq. (|19p in 
the effective theory is 



O = ^({a^nKf^n*. 



(20) 



which is gauge invariant under collinear gauge transfor- 
mations. The sum of all graphs in the full theory with n\- 
collinear gauge emission off any of the particles 1, . . . , r 
in the full theory operator Eq. (p~9|) is equivalent to n\- 
collinear emission from or from the Wilson line 

12| in the operator Eq. (f2T)]) . 
The structure Eq. (|20p is non-trivial, and requires com- 
bining terms with gluon emission from all the particles, 
and using the fact that the operator is a gauge singlet. 
The Feynman rules for multiple gauge emission of n\- 
collinear gluons from particle i gives factors of the form 



£ ■ rij 
k ■ rij 



(21) 



The ni-collinear gauge field has momentum k and po- 
larization e in the ni-direction at leading order in SCET 
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4 3 4 3 4 3 




Pi P2 

FIG. 4: Graphical representation of ni-collinear interactions 
in SCET. 



power counting, so the above expression can be replaced 
by 



FIG. 5: Pair production q(pi) 
runs vertically. 



+ q(p2) -> q'(p4)+q'(p3)- Time 



e • ni 



k ■ rii 



ni ■ tii k ■ n\ k ■ n\ 



(22) 



using the leading (first) term in Eq. ([3]) for the decom- 
position of both k and e. This expression is independent 
of rii . This means that one can change the direction m , 
provided rii ■ n\ remains leading order in the power count- 
ing, i.e. rii does not become almost parallel to n\. One 
can thus move all the rii labels so that they all point in 
a common direction, which can conveniently be chosen 
to be fi\. This choice only makes reference to particle 1, 
and has no information about the directions of the other 
particles. This is the basis for soft-collinear factorization. 

In this paper, we will use the analytic regulator (2f| 
I2H ] used in CGKMl, CGKM2. With analytic regularization, 
Eq. (J22J) becomes 



e ■ n\ 



(k ■ rii) 



1+6 



, -,1+5 n - -,1+S 

(ni ■ rii) \k-m) 



e ■ ni 



(ni ■ ni) S (k ■ ni) 1+S 



(23) 



and the n$ dependence no longer cancels. Thus the iden- 
tities which allowed one to combine all the ni-collincar 
emissions into a single Wilson line in the fi\ direction 
no longer hold. This is a big drawback of the analytic 
regulator. It is possible to use other regulators which 
do not have this problem (35|, but then there are other 
subtleties which must be addressed, related to zero-bin 



subtractions 
factorization 




which are necessary for soft-collinear 
IHH. With the analytic regu- 
lator, Tii-colllinear interactions cannot be encoded in a 
single Wilson line in the fi\ direction; instead one needs 
to include Wilson lines along the directions of all the 
other particles. In the scattering case, this means that 
ni-collinear interactions at one loop are given graphically 
by Fig. 2J This is equivalent to evaluating the collincar 
graphs in the full theory using the method of regions with 
an analytic regulator. We have followed this procedure 
because it allows for a direct comparison of our interme- 
diate results with previous work. 



V. SUDAKOV CORRECTIONS TO 
SCATTERING PROCESSES 

In this section we use the toy model to calculate the 
amplitudes for qq — > qq, qq — > qq, qq — > ti, and qq — > qq c , 
where q denotes a colored scalar particle such as a squark. 
We will call the gauge symmetry color and the particles 
quarks. The corresponding results in the standard model 
are given in Sec. IVIII 

An interesting result is that the SCET 5-matrix el- 
ements are given by summing the results for the two- 
particle case, the on-shell Sudakov form-factor given in 
CGKMl, CGKM2 over all pairs of particles. We first com- 
pute the qq — > qq amplitude explicitly by summing the 
diagrams, and show how the answer can be written as 
a sum over two-particle S'-matrix elements. The general 
proof is given in Sec. IV El 

In this section, as in CGKM2 we use the decomposition 



G = CW + —C^ 

47T 



(24) 



of coefficients and anomalous dimensions into their tree- 
level and one-loop values. In the next section on the stan- 
dard model, we will explicitly include the a /(An) factor 
in the defintion of since there are several different 
gauge coupling constants. 



A. Light Quark Production 

We start with light quark pair-production, qq — > q'q' . 
The kinematics for qq — > q'q' is illustrated schematically 
in Fig. where the incoming and outgoing particles have 
momenta pi, P2 and p^, p±, respectively, and we work in 
the limit s,t,u^> M 2 3> mf. The external particles are 
all on-shell (pf = mf). The Mandelstam variables are 
s = (Pi + P2) 2 , t = (p 4 - pxf and u = (p 3 - pi) 2 . We 
assume q and q' are different flavors, so that only the s- 
channcl annihilation graphs contribute. Identical flavors 
are discussed in Appendix [Al 

At the scale fj, ~ Q the full theory is matched onto 
SCET, and the full theory amplitude at leading order in 
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the power counting is expressed as a sum of local opera- 
tor matrix elements, as in Eq. |T]). The gauge- invariant 
operators in the effective theory are 



0, L L 

Ft R 



Or. L L 

Ft n. 



(25) 



There are only two operators which contribute because 
the fermions are in the fundamental representation of 
the gauge group. For other representations, there can 
be more invariants which contribute, e.g., for isospin one 
fermions, there are three invariant amplitudes in the I = 
0, 1, 2 channels. 
At tree-level, 



MO) 



C 



(0) 
1LR 



r 



(o) 

1RL 



(o) _ 4?ra 

^1RR — 



c. 



(o) 

ILL 



r 



(0) 
2LR 



C i 



(0) 
2RL 



c, 



(o) 

2RR 







(26) 



from the graph in Fig. fl}. 

The one-loop corrections in the full theory are given 
by the diagrams in Fig. [6j as well as vacuum polarization 
and wavefunction graphs. The one-loop corrections in the 
effective theory are given by computing radiative correc- 
tions to the matrix elements of the four-fermi operators 
Oi with tree-level coefficients, and the one-loop match- 
ing corrections C- 1 ' are given by the difference of the two 
computations. The graphs in the effective theory vanish 
on-shell in dimensional regularization, so the one-loop 
matching coefficients a re g iven by the full theory graphs 
computed on-shell [13, l4ll. |42| . In the full theory match- 
ing computation, infrared scales such as the gauge bo- 
son mass M and fermion masses m,, which are all much 
smaller than Q, can be set to zero. Thus the coefficients 
Ci are given by the graphs in Fig. ^ with all masses set 
to zero. The computation is summarized in Ap pendix [Al 
and agrees with previous calculations 0, l43l |44| . The 
one loop coefficients are (removing an overall a/(4ir), see 
Eq. (El): 



C 



(i) 

ILL 



C 



(1) 
1RR 



Ana 



X i s i t ) a AM) 



a 



(i) 

2LL 



C 2RR = --^-/(M)Cl 



C 



(1) 
1LR 



(1) _ 4?ra 



X(s,u) 



f(s,u) 



a 



(i) 

2LR 



c. 



(1) 

2R.L 



,s 



(27) 



where 
X(s,t) 



= 2C F -L 



3L S + — - 8 





FIG. 6: One loop corrections to pair production in the full 
theory. Wavefunction and vacuum polarization graphs are 
not shown 



+C A \^L 2 S - 2L 



-s— t^-s L s 



— L s I Tpup 

3 9 ' 



85\ 
~9~ 



9 I T * n > 



f(s,t) = - 



2s s(s + 2t) , l2 

4L s L t /(_ s _ t ) . 



(28) 



Here rip and ns are the number of Dirac fermions and 
complex scalars. The group theory invariants Cd and C\ 
are defined in Eq. (|40|) , (j4Tj) below. The high scale match- 
ing is the only piece of the computation which cannot 
be obtained by summing the Sudakov form-factor results 
over all pairs of particles. 

If the initial and final quark flavors are identical, then 
there are also ^-channel graphs which contribute to the 
matching (see App. |A")) . 

The next step is to compute the anomalous dimension 
in SCET between Q and M, and the matching corrections 
in SCET at M when the gauge bosons are integrated out. 
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Both results can be obtained simultaneously by comput- 
ing the on-shell matrix elements of Oi in SCET. The fi- 
nite part of the graph gives the matching correction, and 
the infinite part gives the anomalous dimension. The 
SCET diagrams are rii-collinear diagrams and ultrasoft 
graphs. As in CGKMl, CGKM2 the ultrasoft graphs vanish 
on-shell with the analytic regulator, so the only graphs 
which contribute are the collinear graphs. 

The one- loop n^-sector graphs are given in Fig. |4|. 
Particle i is given by the field , and the remaining par- 
ticles are represented by Wilson lines. The computations 
are done using the same analytic regularization method 
used in Refs. CGKMl, CGKM2. The regulated nj-collinear 
propagator denominator is 



1 



2\5i 



2ll+<5. ' 



(29) 



The propagator denominator for particle j interacting 
with ni-collinear gluons becomes 



1 



l(Pj + k)'- 



[2 Pj ■ k] 



(30) 



At leading order in SCET power counting, pi and k are 
ni-collinear, so = nf (fii ■ pi) /2, fc M = (hi ■ k) /2 
and 



1 



( Pl +k) 2 [Km ■ Pj )(n j -n i )} 



l+Si 



(31) 



Thus the analytic continuation of the Wilson line propa- 
gator arising from particle j is 



1 



rii ■ k 
,0') 



(ni ■ ky+^ 



\{n 2 -Pj)(nj ■ rii)] 



(32) 



The key observation is that the Vj regulator parameter 
when particle j is the rij-collinear field £, nj ,pj is related 
a) 

to the v\ regulator parameter when particle j interacts 
with rij-collinear gluons as a Wilson line. This feature 
was already studied in CGKMl, CGKM2 and leads to a 
calculable logarithmic violation of factorization, as dis- 
cussed further in Sec. I VII 

The rii collinear graph with the particle j Wilson line is 
then identical to the n\ collinear graph interacting with 
the 7i2 Wilson line result in CGKMl, CGKM2 with the re- 
placement v\ — > Vi for the collinear particle regulator, 
and i*2 — > v\ for the Wilson line regulator. The regula- 

tor variables Vi,v\ only appear in logarithms, and v\ 
only appears in the boost-invariant combination 



rii ■ Pi 



r (rij ■ Pj )(rij - rii) (m ■ Pi ) 2pi ■ Pj 



(33) 





FIG. 7: Collinear graphs involving particles i and j which are 
related. 



In the Sudakov form-factor results in CGKMl, CGKM2 2pi ■ 
Pi = Q 2 , and Eq. (f3"3")) was the origin of the logQ 2 terms 
in SCET. Here 2pi-pj depends on the kinematic variables, 
and gives a dependence on logs, logt and log it. 

In the Sudakov form-factor computation, there was a 
non-trivial cancellation between the n-collinear and n- 
collinear graphs, so that the sum of the graphs was inde- 
pendent of the analytic regulator parameters i>\. There is 
a similar cancellation here. There are two graphs which 
are related to each other: graphs with gauge boson ex- 
change between i and j in which i is ni-collinear and j 
is a Wilson line, and in which i is a Wilson line and j 
is rtj-collinear (see Fig. [7]). These graphs have identi- 
cal color factors. The regulator cancellation depended 
on two identities given in Appendix A in CGKM2. The 
corresponding relations here are 

( v U) » 2 \ ( ^\ 

log r- 5 — + log -V - log rr 1 — + log 4 



rii ■ Pi M 
2v\v' 1 

log ~ ~ — \" 

/i 2 (rij - pj) (rii - Pi) (rii ■ rij) 



Pj M 




= log- 



2p z v 



2,, 2 



v i i n i -Pj)( n i - Pi) (ni - rij) 
= 2 log 



+ (i <-> j) 



= 2 log 



(rij - pj) (ni-pi)\(m - rij) 
M 2 



2pi ■ pj 



(34) 



which follow from Eq. (f3"2"| . so the v cancellation con- 
tinues to hold. Thus the collinear graphs arc obtained 
by the collinear graphs in the Sudakov form-factor case 
with the replacement Q 2 —> 2pi ■ pj , and summing over 
pairs with the appropriate group theory factor. The ul- 
trasoft graphs vanish on-shell, as in the Sudakov form- 
factor case, so the complete answer is given by adding 
the wavefunction renormalization graphs to the collinear 
contribution. 

The terms which depend on log(pi • pj) arise from the 
regularization of Wilson lines using the analytic regula- 
tor. They depend on the momenta of both particles, so 
it is clear that in ni-collinear graphs, it is not possible 
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to combine the Wilson lines for the other particles into a 
single Wilson line, as that would lose information on the 
Pj dependence. 

Note that the r-particle result obtained by combining 
the Sudakov form-factors over all pairs of particles is valid 
even if all the momenta flowing into the operator do not 
add to zero, i.e. even if there is some momentum inserted 
at the vertex. In the case of 2-particle scattering, we are 
interested in operator insertions at zero momentum, and 
the six pi ■ pj invariants can be written in terms of two 
independent Mandelstam variables. 

The SCET graphs do not change the Lorentz or chiral 
structure of the operators, and only cause rearragements 
of the gauge indices. Thus Gill can mix only with 2 ll- 
Furthermore, the mixing matrix for Oihh' , 2 hh' is inde- 
pendent of the chirality labels h, h'. To keep track of the 
gauge indices, it is convenient to denote 0\hh>, 2 hh' by 



O x = t a ®t a 

e>2 = 



(35) 



The SCET graphs are then a 2 x 2 matrix in 0\^ 2 space, 
and a unit matrix in chirality (h, h!) space. 

The sum of the n-collinear and n-collinear vertx graphs 
with the gauge factor Cf omitted is 



r(Q 2 ) = 



Q 

47T 



-4L 



M 



e e 
t 5ir 2 



-M 



2L M L 



M<-Q 



(36) 



The wavefunction renormalization, omitting group the- 
ory factors is 5 



5Z = T 

47T 



(37) 



The sum of graphs in Fig. (0) which connect particles 1 
and 2 is thus 



r 12 (-2p rP2 )fW®f 

if the operator at the vertex is 0\ , and 
r 12 (-2pi -p 2 )t b t h ® 1 



(38) 



(39) 



if the operator at the vertex is 2 . The minus signs 
relative to Eq. Q36p arise because both momenta pi j2 are 
incoming, whereas in Eq. (|36[) computed in CGKM2, p\ 
was incoming, p 2 was outgoing, and Q 2 — 2p\ ■ p 2 . It 
is useful to add subscripts to T denoting the particles 
involved in the diagram. 



The group theory factors can be simplified using 
t a t a = C F 1 



t a t b t a 



Cf — ~jpA S '"' 



t a t° ® t a t° = d 1 ® 1 + - (C d - C A ) t a <g) t a 
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t a t°®t°t a = Ci \ ®l + -{C d + C A )t a ®t a (40) 
in the notation of Ref. |4^|. For an SU(N) gauge theory, 
C A 
C f 

Ci 



N 




N 2 - 


1 


2N 




N 2 - 


4 


N 




N 2 ~ 


1 



AN 2 



(41) 



so C F = 4/3, C A = 3, C d = 5/3, C x = 2/9 for SU(3) 
and C F = 3/4, C A = 2,C d = 0, C x = 3/16 for 577(2). 
The matrix element of 0\ is 

C f -\ca^ (t a ®t a )T 12 {~2 Pl -p 2 ) 
f-\c, 
L l® 1 + - 
L l® 1 + - 
l 1® 1 + - 
L l ® 1 + - 



(42) 

The terms are given by summing over the six possible 
choices of particle pairs, and including the wavefunction 
contribution for each particle. The terms from gluon ex- 
change between 13 or 24 have minus signs, from charge 
conjugation, since both lines have color flowing into the 
vertex. 6 

For 02, one has instead 

C F (1 <8> 1) r 12 (-2 Pl • pa) + C F (1 <8> 1) r 34 (2p 3 • Pa) 
+ (t a ® t a ) r 14 (-2 Pl • Pi ) + {t a ® t a ) r 23 (2p 2 • pa) 
- (t a ® t a ) r 13 (2pj • pa) - {t a ® t a ) r 24 (2p 2 • Pi ) 

-- (SZi 1 + SZ^ 1 + SZ^ 1 + 5Z1 1 ) C F (1 ® 1) . (43) 



Pi) 




r u (2pi 


■Pi) 


r 23 (2p 2 


■Pa) 


ri 3 (2pi 


■Pa) 


r 24 (2p 2 


■Pi) 


C F {t a <S 


t a ) 



5 These are Eq. (43) with the wavefunction correction removed 
and Eq. (40) of CGKM2. Since we work on-shell, there is no 
need to introduce infrared modes whose virtuality is governed 
by the off-shellness |45H . 



6 Equation [142 t i is true even if there is non-zero momentum in- 
serted at the operator vertex, so that pi + P2 ^ P3 + P4- 
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Equations (|42|) . (|43|) can be written in matrix form, by 
defining the matrix 



tz 
n 



HI + lis 

c F 



r 12 (-2 Pl • pa) - idZf 1 -^ 2 _1 



+c F 



r 34 (-2p3-P4)-^^ 1 -^r 1 



ft* 



ri 



^2 = 



\C d r 1 + \C A r2 ri 



Cm 







r 14 (2pi • p 4 ) + r 23 (2p 2 • p 3 ) 

-Fi3(2pi • p 3 ) - r 24 (2p 2 • Pa) 

ri 4 (2pi -p 4 ) + r 23 (2p2 -p 3 ) 

+ri 3 (2pi • p 3 ) + r 24 (2p 2 • Pa) 

-2r 12 (-2 Pl • pa) - 2r 34 (-2p 3 • p 4 ) . (44) 



Equation (j4"4"|) has an interesting structure — It has a di- 
agonal piece TZ, which is the sum of the on-shell Sudakov 
form factor graphs (including wavefunction factors) for 
1 — > 2 and 3^4, and a term TZs, which depends on the 
amplitude linear combinations r\ and r 2 . TZs contains 
differences of T^- . One can include wavefunction factors 
in TZs by the replacement 

r 

1 1.1 



Sij = r« - \hZ? - X -bZ^ (45) 



without changing t\ and r 2 . We will thus use Eq. (|44|) in 
the form 



TZ 
TZ 

n s 

n 
ri 



TZt + TZs 

CfSi 2 {-2pi -p 2 ) + Cf5 34 (-2p 3 -p 4 ) 
\C d r x + \C s r 2 n 

Cm 



Si4(2pi • pa) + S 23 (2p 2 ' P3) 
~5i 3 (2pi • p 3 ) - 5 24 (2p 2 • p 4 ) 

5i 4 (2pi • p^ + S 23 (2p 2 • P3) 
+5i 3 (2pi • p 3 ) + 5 24 (2p 2 • p 4 ) 
-2«S 12 (-2p x .pa) - 25 34 (-2p 3 -p 4 ) 



(46) 



where 5 is the on-shell Sudakov form-factor including 
wavefunction corrections, i.e. an S'-matrix element, with- 
out any color factors. The r% and r 2 terms contain dif- 
ferences of Sudakov form- factors, and so do not contain 
Sudakov double-logs, which are universal, don't depend 
on particle type, and cancel in the difference. 

The on-shell matrix element of the effective Lagrangian 
CiOi including wavefunction factors is 



<Oi) (0) (O 2 ) (0) 



(i + n) 



Ci 

c 2 



(47) 



where Ci are the operator coefficients and are the 
tree-level matrix elements. 

Equations (j46]l . (|47|) are master equations we will use 
for our scattering computations. For example, to com- 
pute the matching correction when the massive gauge 
bosons are integrated out, we use 



Ci 

c 2 



= (1+fl) 



C 2 



(48) 



where R is the finite part of TZ, and C and C are the 
coefficients in the high-energy theory with gauge bosons 
and the low-energy theory without gauge bosons, respec- 
tively. Similarly, the anomalous dimension matrix is 



"Cr " 






.C2_ 


= 7 


.C2 



(49) 



where 7 is the anomalous dimension computed using the 
l/e terms in TZ, i.e. —2 times the 1/e terms in TZ at 
one loop. The matching conditions and anomalous di- 
mensions are given by Eq. (|46|) with Sij replaced by the 
corresponding Sudakov form-factor matching correction 
and anomalous dimension computed in CGKM2 without 
any additional Fcynman graph computations. 

We now apply the master formula to the SCET anoma- 
lous dimension for qq — > q'q' in the region Q > /1 > M, 
and to the matching condition at M. The anomalous di- 
mension is given using Eq. (|46p with S replaced by the 
SCET anomalous dimension for the Sudakov form fac- 
tor, i.e. by for the bifermion operators in Table I of 
CGKM2, S — > 4Ln — 6. The anomalous dimension matrix 



is 



7?> 



2C F ^4 log— -6 

~ 2C d \og± + 2C A \ogf 
8Cilog£ 



'log- 




(50) 



or using the notation defined in Eq. (U), 

7 (D = ^(i) 1 + 7 W 
= 2C F (4U-6) 



7 V 



7* W 



2C d l t 



2CA\-ut/s 2 8L t / u 



8C1 L t / M 







(51) 



All logarithms of negative argument are defined by the 
branch log(— s — i0 + ), and \og(ut/s 2 ) = log(— u — i0 + ) + 
log( — u — i0 + ) — 21og(— s — i0 + ), etc. as discussed ear- 
lier. The off-diagonal terms vanish at t = u, i.e. when 
the center-of-mass scattering angle is tt/2. 7s is called 
the soft anomalous dimension. We will see explicitly that 
the soft anomalous dimension and the soft-matching Rs 
are universal, and independent of the external states, i.e. 
they are the same for fermions and scalars, and indepen- 
dent of the particle masses. In our computation using 
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the analytic regulator, 75 arises from collinear graphs; 
the ultrasoft graphs all vanish. 

The anomalous dimension 7 is twice the anomalous 
dimension for the Sudakov form-factor, and contains a 
log(— s//j, 2 ) term which produces double logs in the am- 
plitude on integration. The soft anomalous dimension 7s 
does not contain any parameterically large logarithms, 
since s,t, u are all formally of order Q 2 . 

The matching matrix is given by replacing S by the 
matching for bifermion operators in Table I of 

CGKM2 S 



-L 2 M + 2l M l Q - 3L M + 9/2 - 57H/6, 



2C F 



-L 2 M + 2Lj\/L s — 3Lm + - 



G 



R 



(i) 



-M 



CdS-t/u + CA^-ut/s 2 4L t /„ 



4CiL 



t/% 



(52) 



Note that there is a non-trivial low-scale matching cor- 
rection. At [i — M, Lm — 0, and Rg vanishes. This is 
an accident in the toy model at one-loop. In the standard 
model, R^ does not vanish, and has terms of the form 
log /Ml . 

This completes the computation for quark production. 
The matching and anomalous dimensions are combined 
to give the final amplitude in the usual way, and give the 
exponentiated SCET form for the Sudakov logarithms 
discussed in CGKM2. The matrices 7s and Rs are univer- 
sal, and have the same values for heavy quark production 
and for squark production, as we see explicitly below. 



B. Light Quark Scattering 

The next process we consider is light quark produc- 
tion, q(pi) + q(p 3 ) -> q'(p 2 ) + q'{pi) (see Fig. [5J with 
q =/= q' , which is related to quark scattering in Fig. [5] 
by crossing symmetry, with the replacements p 2 — + — P2, 
pj, — > —p3- The Mandelstam variables for quark scatter- 
ing are s = (pi + p 3 ) 2 , t = (p 2 - Pi) 2 and u = (p 4 -pi) 2 , 
so the amplitudes are obtained from those in the previ- 
ous section by the replacement s — > t, t —> u, u — > s. 
Identical flavors are discussed in Sec. [5] 

The anomalous dimension matrix is 



7 



(i) 



7?> 



= 7«1 
= 2C F (4L t -6) 

2CdL„/ s + 2CA\- us /t 2 
8CiL u / s 



8L 



u/s 



(53) 



and the matching matrix is 



R (1) = 2C F 



R 



(1) 



2L M L t -3L 



5^ 
"6" 




FIG. 8: Quark scattering q(pi) + q'ips) — > qijfr) + q'(pi) 
Time runs vertically. 



R 



(i) _ 



-M 



Cd)-u/s + CA^-us/t 2 4L M / S 

4CiL u / s 



(54) 



C. Heavy Quark Production 

Consider the annihilation of a light quark-antiquark 
pair to produce a heavy quark-antiquark pair, sugges- 
tively labeled tt, via the process q(pi) + q{p2) — * t{pi) + 
i{ps). The kinematics and Mandelstam variables are the 
same as Sec. IV Bl the only difference is that the final 
particles have mass m which is not negligible compared 
with the gauge boson mass M, but is much smaller than 
Q, so that s,t, u - Q 2 > m 2 , M 2 . 

The first step is to match the full theory onto SCET 
at pi ~ Q. The fields £„ 3 and ^ n4 are now taken to have 
mass m [13, EH . The matching condition at Q can be 
computed by from the full theory graphs with all scales 
much smaller than Q set to zero, so the matching at Q 
is the same as for the light quark case. 

The second step is to run SCET operators in the effec- 
tive theory from Q to m. The SCET anomalous dimen- 
sion is independent of low mass scales and again gives 
the same result as in the massless case, Eq. (j5"3"]) . 

The third step is to switch at the scale fi ~ m to an 
effective theory where the heavy quarks are described 
by heavy quark effective theory (HQET) fields t V3 and 
t V4 [H. The four-fermi SCET operators of Eq. are 
matched onto the SCET/HQET operators: 

0i ->0J - i v J a ^P L J V3 [C n2 W n2 }t a ^P L jW^ ni } 

2 ^0' 2 = t Vi ~f»Pit V3 [£ n2 W n2 ] lfi PL [Wltm] ■ (55) 

The HQET fields do not transform under a collinear 
gauge transformation; therefore, there is no factor analo- 
gous to the Wilson line that goes along with £ n . The 
heavy fields t Vi still couple to ultrasoft gauge bosons. 

The matching condition at /1 ~ m is given by com- 
puting the difference of the graphs in the theory where 
particles 3 and 4 are described by SCET fields, and the 
same graphs computed when the two particles are de- 
scribed by HQET fields. Particles 1 and 2 continue to be 
described by SCET fields. The group theory and kine- 
matic factors for each pair of particles remain unchanged 
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as we switch from SCET to HQET, so the matching con- 
dition is given by Eqs. (|4*6")) with each T being replaced 
by the difference of the corresponding graph in the two 
theories. Thus one can use 



Sia -> 
S34 — > Rhh 
Sij{x) -> Rm 



ij = 13, 14, 23, 24 



(56) 



where R^h is the matrix element for the Sudakov form- 
factor in going from two SCET to two HQET fields, and 
Rhi is the matrix element for the transition from two 
SCET fields to one SCET and one HQET field, dropping 
any overall group theory factors. The matching coef- 
ficients can be read off from Eq. (80) and Eq. (85) in 
CGKM2 



S\2 


-> 






-> R 


+ T 




-> R 


1 


R 


= T 



Li 



13,14,23,24 



1. 



7T" 

12 



(57) 



where we use the entries from the first rows of Tables II 
and IV of CGKM2. Thus 



K 



(i) 



n w 1 + 

C F (R + T) 
C F [ Li - L„ 





7T 

IT 



(58) 



using Eq. (|46|) for the matching. vanishes since n = 
2R-2R = and r 2 = 4R - 2{R + T) = 0. 

The anomalous dimension below m is given by using 
Eq. (|46l) with the replacement analogous to Eq. |57|) for 
the anomalous dimension, 



5i 2 (-2pi 


P2) - 


■* 7i (- 


-a) 


5 3 4(-2p 3 




* 73 (- 


s) 


5i 3 (2pi 


P3) ~ 


72 (- 


-u) 


5i 4 (2pi 


Pi) ~ 


72 (- 


-t) 


5 2 3(2p 2 


Pa) ~ 


72 (" 


-t) 


S 2 i{2p 2 


Pi) ' 


■+ 72 (- 


-u) 



(59) 

where 71,2,3 (Q ) are the entries from the first rows of 
Tables I, II and IV of CGKM2. They are the anomalous 
dimensions for 11, hi and hh currents, respectively. The 
anomalous dimension matrix in the HQET/SCET theory 
is 



7 



(i) = 

? W = 

ri -> 

f2 ~> 



7i 1} 



7 V *' 1 
C F ( 7l (-s)+7 3 (-s)) 
2 72 (-i) - 2 l2 {-u) 

2 72 (-t) + 2 72 (-m) - 2 7l (-s) - 2 73 (-s) 



71 (Q 2 ) 

72 (Q 2 ) 

73 (Q 2 ) 

r(u>) 



4L Q -6 

4Lq - 2L m - 5 

4 [wr(to) — 1] 



log (u? + yjw 2 — l) 



1 



2m? 



1 



(60) 



Since we are working in the limit Q 2 ^> m 2 , wr(w) — 1 — > 
log(2iu) — 1 — > log(<5 2 /m 2 ) — 1 up to power corrections. 
This gives 



7 V 



7W1 



7^ = 



C F (8L s - 4L m - 10) 
2C , dL t / u + 2CUL ut / s 
8CiL t / u 



8L 







(61) 



The last step is to integrate out the gauge boson at fi ~ 
M and transition to the theory with no gauge bosons. 
The matching is given by Eq. (146|) where Sij are replaced 
by the corresponding results for the Sudakov form-factor 
matching, 



S 12 {-2p\ -p 2 ) 
S 34 (-2p 3 -p 4 ) 
S 13 (2p 1 -p 3 ) 
S 14 (2pi -p A ) 
S 23 {2p 2 -p 3 ) 
S 2i (2p 2 - Pi) 



D(-s) 
U(-s) 
S(-tt) 
S(-t) 

s(-t) 

S(-u) 



(62) 



where D, S, U are given in the first rows of Tables I, II 
and IV, respectively, of CGKM2. The matching is 



rW = R^t+R 
R^ 

ri - 
r 2 

so that 

R^ 1 

RW = Cf 



C F (D(-s) + U(-s)) 
2S{-t) - 2S(-u) 

2S(-t) + 2S(-u) - 2D(-s) - 2U{-s) (63) 



= R {1) t + R 



(i) 



M 



-5L M + - 



4L M L 

6 



2Lj\/L„ 



pC 1 ) _ , 



M 



Cd\-t/u + Ca\-ui/s 2 4L t / u 
ICxU/u 



(64) 



In summary, the computation proceeds as follows: (a) 
Match at fj, - ^/s using Eq. J25]), (HH) (b) Run be- 
tween s/s and to using Eq. (|50|) (c) Matching at to using 
Eq. (fS"8")) (d) Run between to and M using Eq. (|6"T|) (e) 
Match at M using Eq. 
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(a) 



(b) 



FIG. 9: Tree level squark production in (a) the full theory 
and (b) the effective theory. 



If the fermion mass is not much larger than M, as 
is the case for the top-quark, one can replace (c), (d) 
and (e) by a single step, (c') Integrate out the fermion 
and gauge bosons simultaneously at fi ~ m ~ M, as in 
Sees. VIII D,G of CGKM2. In this case, the matching is 
given by Eq. (87), (91) of CGKM2: 



5i2(-2p!-p 2 ) -> D(-s) 

S M (-2p 3 - Pi ) -» D(-s)+2f F (z)-h F {z) 

S l3 (2 Pl -p 3 ) -» D(-u) + / F (z)-ft F (^/2 

5i4(2pi-p 4 ) -> £>(-*)+ /p(«)-/»f(«)/2 

5 23 (2P2-P 3 ) -» D(-t)+f F {z)-h F {z)/2 

S 2 i{2 P 2-Vi) -> D(-u) +/*.(«)- M*)/2 

2 

m 



(65) 



where the functions /i? and hp are given in Appendix B 
of CGKM2. They are the change in the matching condition 
due to the quark mass. The matching matrix becomes 



RW t + R 



(i) 



Cp(2D(-s) + 2f F (z)-hp{z)) 
2C F ( 

+ f F (z)-h F (z)/2) 

Cd\-t/u + CUL„ t / 
4dL t/u 

The f F and hp terms cancel in Rs- 



R 



(i) 



5^_ 2 
6 



-M 



2LmL s — 3L M + - 



4L t/u 




(66) 



D. Squark Pair Production 

As the final example, we consider heavy scalar (squark) 
pair production via q(pi) + ?(p 2 ) — * ^(p 3 ) + t*{Pi), where 
t, i* are the squark and anti-squark. The squarks are 
taken to have mass m <C This example shows 



how one can use the scalar and scalar/fermion results 
in CGKM2 to compute squark production. The discussion 
parallels that for heavy quark production in the previ- 
ous section. The only difference is that since some of the 
particles are scalars, we need to use the 'tpip, </yif> and $4> 
entries from the tables given in CGKM2. 

The first step is to match onto SCET at the scale fi ~ 
y/s. The four-particle operators are 



02 = [^ ni w ni wv 3 + i v 4 r[wl^ n3 ] 



(67) 



where %D 3 = V + g(n 3 ■ A„ 39 )(n 3 /2), iX> 4 = V + #(n 4 • 
A n4 q)(n4/2) are the label covariant derivatives on parti- 
cles 3 and 4, respectively. 
The tree-level coefficents are 



(0) 



Aira/s 



cf > = 



(68) 



from the graph in Fig. ©. 

The anomalous dimension in SCET below the scale Q 
is given by using Eq. (146|) . and the values for the graphs 
in the region between Q and m given in CGKM2. For 
the anomalous dimension matrix, this means the replace- 
ments 



S 12 (-2pi -p 2 ) 

S 3 4(-2p 3 -p 4 

Si 3 (2p 1 -p 3 ) 
Su{2pi - Pi) 
S 23 (2p 2 -p 3 ) 
5 24 (2p 2 - Pa) 



7w(-s) 
7i^(- M ) 

7iV0( _t ) 
7i#(-«) ■ 



(69) 



The anomalous dimensions are given in Table I of CGKM2. 
The subscript iptp, ipej) and (jxj) means we use the anoma- 
lous dimension for bi-fermion operators, fermion-scalar, 
and bi-scalar operators respectively. The anomalous di- 
mension is 



7 



(i) 

'2 



7« 



7^' a 

Cp (ji^(-s) + 7i00(-s)) 
2Cp (4L S - 7) 

2j2if,4,{-t) - 2"/2ii>4> (-u) = 81o S^ 

272^(-*) + 27 2 ^0(-w) - 27i^(-s) - 27 W (-s) 
ut 

81og^ 



2C,jL t /„ + 2C y iL„ t / s 2 8L t / u 



8C1L 



t/v 







(70) 



After running the operators down to fj, ~ m using 
Eq. ([70)1 , one matches to an effective theory in which the 
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scalars are replaced by HQET fields. This is given by us- 
ing Eq. ([56")) , where the scalar values of R^h and Rhi are 
used. This means in Eq. ([57f , R + T should be replaced 
by the bi-scalar value on the second rows of Tables III 
and IV, and for Ty , R should be replaced by R+t^, the 
entry on the fourth row of Table III for a bilinear with a 
heavy scalar and massless fermion: 



n 

r-2 

(i) 



K {1) 1 + 

Cf (R<j><p + T^) 



C F 




2U 



Tl- 



K 



0. 



(71) 



The running in the HQET/SCET theory below to, and 
the matching at M is identical to Eq. (pUj) and Eq. ([6"3")) 
in the previous section, since it does not matter whether 
the HQET field is a scalar or a fermion. 

In summary, the computation proceeds as follows: (a) 
Match at \i ~ y/s using Eq. (|57|) (b) Run between 
and m using Eq. ([50)) (c) Matching at to using Eq. (|TTj) 
(d) Run between m and M using Eq. (|60[) (e) Match at 
M using Eq. ([63]). 

If the squark mass is not much larger than M, one can 
replace (c), (d) and (e) by a single step, (c') Integrate 
out the squark and gauge bosons simultaneously at (i ~ 
to ~ M, as in Sec. VIII G,D of CGKM2. In this case, the 



matching is given by 










S 12 (-2 Pl 


V2) -> 












P4) 




-*)H 


-2/s(«; 


- M*0 


5i3(2pi 


P3) 




-«)- 




-h s {z)/2 


5i 4 (2pi 


P4) -> 




-*)H 


- M*) 


-h s (z)/2 


*5 2 3(2p 2 


P3) -> 




-*)H 




-h s (z)/2 




P4) 
z = 


9 

TO 

M2 


-«)- 


f /s(*) 


- h s (z)/2 



(72) 



where the functions /g and h$ are given in Appendix B 
of CGKM2. The matching matrix becomes 



R.W 



= R^l+R^ 

= C F {D^(-s) + D^(-s) + 2f s (z) - h s (zj) 

3, 



2C f 



-L 2 m 



2L M L 



:1m + 4 



r2 



2D^{-t) - 2D H {-u) = 4L A/ log- 
2D^(-t) + 2D H {-u) - 4L>^(-s) 
4L M log -j 



R 



(i) 



= L 



M 



4CiL 







(73) 



E. Extension to more particles 

In the previous examples, we saw that the four-particle 
5-matrix elements could be obtained by summing the 
two-particle ^-matrix elements over all pairs of particles. 
This result can be generalized to gauge singlet operators 
with an arbitrary of particles. 

The SCET graphs do not depend on the Lorentz struc- 
ture of the operators, the non-trivial dependence is on the 
gauge structure of the operators. We write the operators 
with all incoming fields. An outgoing particle can be rep- 
resented as an incoming field in the complex conjugate 
representation. The incoming fields are combined into a 
net gauge singlet, and we have 



\nt \Oi = o 



(74) 



where T° acts on the indices Oi associated with field a. 
To make the notation clear: Assume ip and x transform 
in the fundamental and anti-fundamental of SU (N) , and 
O = x l ipi- The action of and T£ on O are: 



T«0 = x i (T a ) i j tpj 



T a 



(75) 



Here (T )^ and (T a ) % j are the representation matri- 
ces in the fundamental and anti-fundamental representa- 
tions, so that 



(r a y j = 

from which it follows that 



(T Q ) 



(T$+T%)0 = 0. 



(76) 



(77) 



The sum of graphs with gauge boson exchange between 
particles a and (3, without any gauge factors, will be 
denoted by T a p{2p a ■ pp), as in the preceeding section. 
The graph is computed with momentum p a incoming, 
and pp outgoing. Treating all particles as incoming for 
both color and momentum flow means that the graph 
including color factors is — T a p(— 2p a ■ pp)T®Tp. The 
minus sign of the argument takes care of the change in 
momentum labeling for /?, and the overall minus sign 
is the charge conjugation minus sign from reversing the 
color flow of p. 

The sum of graphs including gauge factors is then 



^-r Q/3 (-2pa-P /3 )(r a % a o i ) 

(a/3) 



(0) 



(78) 



Co) 



where we sum over all pairs (a/3), and (T^T^Oi 

the tree-level matrix element of the operator after the 
action of the gauge operators. 
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The one loop contribution to the on-shell matrix ele- 
ment is 



where A is the scattering amplitude, and F r (Q 2 ) is the 
Sudakov form factor for particle r at some reference mo- 



(0) 



mentum, e.g. Q 



As is called the soft amplitude 



<«/3> 
1 

~2 



in the literature. With this definition, the soft amplitude 
has the form at one-loop 



(0) 



(79) 



(0) 



including the wavefunction corrections for each external 
leg. This can be rewritten as 



]T ~S a p(-2 Pa ■ pp) (Z£Z20,) W , (85) 

<a/3> 



- E 

<a/3> 



1 _i 1 _l 

T a p(-2p a -p p )--SZ a --jSZp 



E 

( a p) 



{T a a T a p O % ) 



(0) 



The first term in square brackets is the on-shell Sudakov 
form-factor for the two-particle case, including the wave- 
function correction, 



11, 

T a p(-2p a -pp) - -6Z a - -SZ p = S a p(~2p a -pp) 



We can simplify the remaining terms using 



= 



E K 1 



SZ: 



(a/3) 



a 

(82) 



since the one-loop Sudakov form factor for particle a is 

( T a T a n \ (a}Ta T S S aa(Q )■ 

\ o ^ '/ The Sudakov form factor has the form at one-loop (see 
the next section) 

S a p(-2p a -pp) = A\og ~ 2pa 2 Pp +B a + Bp, (86) 

M 

(80^h ere A is a universal coefficient independent of parti- 
cle type proportional to the cusp anomalous dimension 
which is known to be universal [50], plus one-particle 
terms B a which depend on the particle type, but are 
independent of p a - pp. 

Using Eq. (Hi with 5Z~ X -> B a , and with SZ^ 1 1 
shows that the soft amplitude is given by a sum of the 
cusp part of the Sudkaov form factors, with coefficients 
which add up to zero, i.e. it can be written as differences 
of A-terms. The B terms all cancel. We have seen this ex- 
plicitly in Eq. (|46|) . Thus the soft amplitude is universal, 
proportional to the cusp anomalous dimension, and for- 
mally has no large log terms since the differences of two 
A terms gives a logarithm whose argument is order unity 
in the power counting, e.g. logt/u = log(— t) — log(— u). 
This also implies that the soft anomalous dimension is 
proportional to the cusp anomalous dimension. While 
the above argument is at one-loop, we believe the gen- 
eral structure persists at higher loops. This property has 
been seen explicitl y at two-loops in a very interesting re- 
cent computation [34j . 



(81) 



which follows from Eq. (|74p . and reduces Eq. (f80l) to 



<a/3> 



(0) 



The final answer can be written directly in terms of the 
on-shell two-particle matrix elements, as we found in the 
previous section for the four-particle case. Equation l[55]) 
is valid even without a summation on the gauge index 
a, and this will be useful in breaking up the electroweak 
corrections into the W, Z and 7 contributions. 

It is conventional to take the multiparticle scattering 
amplitude and divide it by the Sudakov form factors, 



.4 



A, 



VI. FACTORIZATION 

There are strong constraints on the form of the scat- 
tering amplitude in SCET. We will discuss these in the 

(83) context of the analytic regulator. The results hold for 
the S-matrix elements, and so are independent of any 
specific regulator. We have obtained the same results us- 
ing a different regulator [35|. We study the case where 
there is only a single amplitude to avoid problems with 
matrix ordering. This is the case, for example for scat- 
tering in a U(l) gauge theory. If there are several gauge 
structures which can contribute, then the amplitude A is 
a matrix, and one has to worry about matrix ordering. 
For example in SU (N) gauge theory, there are two gauge 
invariant four-particle operators, T a (g) T a and 1 (£> 1, so 

(84) A is a 2 x 2 matrix. We briefly comment on the matrix 
ordering problem at the end of this section. 
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The r-particle scattering amplitudes are given by Uj- 
collinear sectors, i — 1, . . . , r and the ultrasoft graphs. 
With the analytic regulator, the on-shell ultrasoft graphs 
vanish, and we only have to consider the collinear sec- 
tors. 7 

The ni-collinear graphs have the form Fig. IU where 
particle i is given by the SCET field, and all the other 
particles are Wilson lines. The on-shell graph depends 
on the particle masses {mk}, the renormalization scale 
/j, and the analytic regulator parameters. The parti- 
cle masses are the masses of any particles given by rj.j- 
collinear SCET fields, such as the gauge boson masses, 
and the mass of particle i. They can also include the 
masses of other particles which couple to particle i. For 
example, in the standard model, a graph with a final 
ni-collinear i-quark can depend on m t and mi,, since 
n-collinear W bosons couple t to b. The analytic reg- 
ulator parameters for an rti-collinear graph are uf from 

the SCET field, and v*f\ j ^ i from the Wilson lines. 

(i) 

Boost invariance requires v\ to occur in the combina- 
tion vf* /(fli - pi) = Vj/(2p. l -pj), as noted in Sec. El The 
analytic regulator parameters only occur in logarithms, 
and we use the abbreviations 

U = log^ 2 , 
Py = log(2pi-pj) (i^j) 
Pa = 

P-y = Pji ■ (87) 

The total rii-collinear amplitude has the form 

ex P A({L, - Py},{m fe },//,{5 fe }) (88) 

A depends on the momenta of all the particles in the 
process through its dependence on 2pi ■ pj in Ly . The 
total connected amplitude exp A is given by the product 
of the different collinear sectors, so that 

r 

A({pk},{ m k},v) = ^A^iLj -P l;j } ,{m k },fi,{5 k }). 

i=l 

(89) 

The amplitudes A; and A begin at order a. The tree-level 
amplitude is the 1 in the expansion of exp A. 

The individual terms A depend on the regulator pa- 
rameters {vk\ and {5k} and are singular as {5k} — > 0, as 
can be seen explicitly in the one-loop results given in I 
and CGKM2. However, the sum A is finite as {5k} — > and 
independent of the analytic regulator parameters {i>k}- It 
can only depend on the particle masses (including inter- 
nal particles), the external momenta, and /i, as written 
in Eq. ([89]). 

The cancellation of the Lj dependence is a powerful 
constraint on the form of the SCET amplitudes. We 



7 With other regulators, the ultrasoft graphs can be non-zero. 



showed in CGKM2 how it implied that the low-scale match- 
ing D when the massive gauge bosons were integrated out 
had to be at most linear in log Q 2 to all orders in pertur- 
bation theory. In Eq. (|89[) . the right hand side depends on 
momenta only through the terms Py, which occur only 
in the combination Lj — P y in A;. The Lj cancellation 
implies that A and A can be at most linear in P y , to 
all orders in perturbation theory. The proof follows from 
a straightforward but tedious application of the princi- 
ple of separation of variables used in partial differential 
equations — if f(x) + g(y) is a constant, and x and y are 
independent variables, then f(x) and g(y) must both be 
constant. 

The linearity of A in Vij implies that the anomalous 
dimension and low-scale matching conditions D are lin- 
ear in Vij, since they are determined by the infinite and 
finite parts of A, respectively. The only multi-particle 
dependence of A is through the Vij dependence in the 
analytic regulator. Since A is linear in Vij, this leads to 
a two-particle dependence, plus one-particle terms, i.e. 
A has the form Eq. (|86[) to all orders. The A term is uni- 
versal; it cannot depend on the properties of the particles 
such as masses, because it is generated from Wilson line 
vertices which are independent of m. The m-dependcnce 
must be in one-particle contributions. 

If A is a matrix, then the analysis becomes more com- 
plicated, but the general features discussed above con- 
tinue to hold. The SCET anomalous dimension still con- 
tains only a single logarithm to all orders in perturbation 
theory [32|, [HI]. The amplitude Eq. (TIT))) is now a ma- 
trix equation, and the anomalous dimension integration 
is path-ordered in fi. The Lq terms in the anomalous 
dimension 7 are proportional to the unit matrix 1, and 
can be pulled out as an overall multiplicative factor that 
commutes with the non-Abelian exponentiation of the 
integral of the rest of 7. 

The high-scale matching need not be a square matrix, 
and one should replace expC(Q) — * c(Q) in Eq. (fl5|) . 
There are no large logs in either C(Q) or c(Q). The low- 
scale matching D is also not a square matrix. It has the 
form 

d HM))e Dl M M » logQ2 / M2 (90) 

where do is a matrix, and D\ is a number, i.e. it is pro- 
portional to the unit matrix. Thus our result that the 
low scale matching has the form expZ), where D has a 
single log to all orders in perturbation theory still holds 
in the matrix case, in the form Eq. ([90)1 . 

The structure of the amplitudes discussed in this and 
the previous section are a very powerful constraint. They 
follow from renormalization invariance of the effective 
theory and the universality of the cusp anomalous di- 
mension. More extensive comments will be given else- 
where [35:1. 
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VII. APPLICATION TO THE STANDARD 
MODEL 

In this section, we apply the methods developed so far 
to compute radiative corrections in the standard model. 
There are several major differences between the toy the- 
ory and the standard model. The standard model is a 
chiral theory and the couplings of the matter fields to the 
gauge fields are more complicated, with matter fields in 
several different representations of the gauge group. The 
gauge group is not simple and, we have to treat several 
different gauge interactions. After electroweak symme- 
try breaking, there is electroweak gauge boson mixing 
between the W3 and B, which gives W and Z bosons 
with different masses, and a massless photon. Finally, 
there are also Higgs exchange corrections proportional to 
the fermion mass, which are relevant for the top quark. 
It is straightforward to obtain the results for the stan- 
dard model, following the same procedure used for the 
toy model. We have already shown in CGKM2 how to 
obtain the Sudakov form-factor for the standard model 
including all these effects. In this section, we use the 
methods demonstrated in the previous section to calcu- 
late the radiative corrections to dijet, dilepton, top quark 
and squark production in the standard model. These cal- 
culation are a non-trival example of the techniques devel- 
oped in the previous sections and in CGKM2. There are 
eighty independent amplitudes we need to compute, not 
including those related by crossing or flavor symmetry. 

The left handed quarks and right handed quarks are in 
different representations of the unbroken gauge group of 
the standard model. The left-handed quark doublets will 
be denoted by Q ? \ where i = u, c, t is a flavor index, the 
right handed charge 2/3 quarks by U R . , the right-handed 
charge —1/3 quarks by D^, the left-handed lepton dou- 
blets by L L and the right-handed lepton singlets by E R . 

Written in terms of SU (2) components, QW is 




(91) 



and ZyW is 



E 



(92) 



where the neutrinos are weak-eigenstates. 

All the lepton and down-type quark masses can be ne- 
glected in our calculation, so we can work in the weak 
eigenstate basis, the CKM matrix V does not enter the 
SCET computation, and generation number is conserved. 
The only place where V enters is in the matrix element 
of SCET operators in the proton state, i.e. in computing 
the cross-section from the amplitude using the parton 
distribution functions, since these are given in terms of 
mass-eigenstate quarks. 



A. Matching at Q 

For scattering processes, we need to consider four- 
particle operators in SCETew generated at the scale Q 
by the graphs in Fig. Q] and Fig. El The SCET fields cor- 
responding to a standard model field will be given by the 
replacement — ► £ r (Q^)> etc. We will drop the n,p 
labels on £, and instead use the subscript r = 1, . . . , 4 to 
denote the particle label in the scattering process. Thus 
£i(tt) described collinear u quarks with momentum p\, 
etc. 

At the scale Q, the effective Lagrangian is the sum of 
terms representing the scattering of the various particles. 
If the initial and final particles are both quark doublets, 
then the Lagrangian is 



C 



QQ 



a 



x 'HQ {i) )W 2 T A t h ^wl^{Q^)\ 

+c QQ12Ji [UQ U) W a iMUQ U) ) 
x 'UQ {i) W2t a ^wtUQ (l) )] L 
+c QQ21 , fi [UQ u) )W4T A ^w^ 3 (Q U} ) 
'UQ {€> W2T A 1 »wit 1 (Q^)\ L 

+c QQ22Ji [UQ (f) )W 4 iMUQ (f) ) 
'UQ {i) )w 2 ^wU 1 {Q^)\ L . 

We will write this in the abbreviated form 

£QQ = C QQ11Jt [T A t a ] L ®[T A t a ] L + C QQ12Jl [t a ]L®[t a ]L 
+C QQ21Jl [T A } L <g> [T A ] L + C QQ22Ji [l} L ® [1] L . 
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(94) 



The flavor quantum numbers are encoded in the sub- 
scripts on C. Recall that T A are the SU(3) generators, 
and t a are the SU(2) generators. The subscript 1 is used 
for T A <g> T A or t a <g)t a , and the subscript 2 for 1 ® 1. 
Similarly, one has the other terms 
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£-UE 
£-DQ 
£-DU 
£-DD 

£dl 
£de 
£lq 
£lu 
Cld 
Cll 

£-EQ 

£eu 

&ED 
£-EL 
£-EE 



CuE,fi[l] R ® [l]fi 

C D Qiji[T A ) R ® [T A ] L + C DQ2Ji [l} R ® [\\ L 
C D mji[T A ] R ® [T A ] fl + C DMi/i [l] ii ® [l] fl 

CDDi,/,[r A ] fl ® [T A ] fl + cw/*[i]i? ® 

Cx)B,/i[l]je® 

Glqi^Il ® [i a ]i + C £Qa ,/i[l]i ® [1]l 
C.ti7,/i[l]t ® 



Ci2j,/i[l]i [1] R 

C LL1 , fi [t a ] L ®[t a ] L 
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Cj3a,/i[l]ii® [1]« 

CfeL,/i[l]fl ® 
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■C LL2 ji[i\ L ® [i]i 
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The tree-level matching coefficients from the graph in 
Fig. [T] are (f ^ i) 
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and the one loop matching coefficients are 
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for N = 2,3 for SU(2) and SZ7(3), respectively. 
n F N (ns N ) denotes the number of Weyl fermions and 
complex scalars in the fundamental representation of 
SU(N). The matching is symmetric between initial and 
final fermions, so that Cjjq,i = Cqu.i, etc. The coeffi- 
cients Cqdj are given by Cquj with Yjj — ► Yd, Cuu,j by 
Cud, 3 with Y D — > Y v , C D D,j by Cr/pj with Y v ~*Y D , 
Cdl.j by C[/lj with Yu —>Yd, &ndCDE,j by Cus.j with 
Yu — > Yd, and so have not been listed above. For identi- 
cal particles (i.e. CQQj—i, CuL,f=i, etc.) there is also the 
crossed-channel contribution as discussed in Appendix [XI 



The above matching coefficients do not include Higgs 
exchange contributions. The Yukawa couplings are pro- 
portional to the fermion masses, and the only Yukawa 
coupling large enough to be relevant is the top quark 
Yukawa coupling. Higgs corrections only arise at one- 
loop for LHC processes, since the initial state is pp, and 
contains no t-quarks. 8 The Higgs contributions to the 
matching for operators containing in the final state 



arc 1 



One can always treat the proton as a hadron in QCD with all 
heavy flavors integrated out. Heavy quark distribution functions 
in the proton are calculable in terms of light-quark distribution 
functions; see e.g. Ref. [5lll . 
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whereas the contribution matching for operators contain- 
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ing !/(*) in the final state are 
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with i = u,c for Q or £/, and i = d, s,b for Z). The 
logarithmic terms for SCqe.u and 8CjjE,ti agree with 
Ref. [13. 

Once we match onto SCET, Higgs vertex corrections 
are power suppressed, as shown in Ref. 0], and the only 
Higgs contributions in SCET are wavefunction renormal- 
ization corrections. 



B. Anomalous Dimension below Q 

The anomalous dimensions in SCET between Q and m 
are obtained using the results of Sec. IV Al The anoma- 
lous dimension due to gluon exchange depends on the 
color quantum numbers of the initial and final fermions. 
If both are color triplets, then the operators have the the 
color structure Ci T A ® T A + C 2 1 ® 1 . The anomalous 
dimension is given by Eq. ([50)1 with group invariants re- 
placed by their values for N = 3, and with the a/ (4%) 
prefactor for QCD, 



lSU(3) 
lSU(3) 

lS,SU(3) 



lSU(3) 1 + 1S,SU(3) 

8 a 3 



3 4tt 
4tt 



(4U - 6) 

ir'-t/u + 6L ut / s 2 8L t / u 
T l t/u 
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This 2x2 anomalous dimension matrix acts on operators 
with color structure T A ®T A and 1(8)1, and does not mix 
different flavors, chiralities or SU(2) quantum numbers. 
Thus the renormalization group equation has the form 
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where (Ci,C2) are the pairs 

(CQQ21.fi, CQQ22Ji), 
(CQDl,fi,CQD2Ji), 
(CuUIJi, CuU2Ji), 
{CDQl,fi,CDQ2Ji), 
(CDDl,fi,CDD2J'i)- 

If one of the fermions is a color triplet and the other is a 
color singlet, the operator has the color structure C 1(3 1. 
The QCD anomalous dimension for C is identical to the 
Sudakov form-factor case, 



(CQQllJi, C QQ i2Ji), 
(CQUIJi, CQU2.fi), 
{CuQl,fi,CuQ2,fi), 
(CuDIJi, CuD2,fi), 

(Cuu\,fi,Cuu2ji), and 
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3 4tt 
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If both fermions are color singlets, then 



1SU(3) 



= 0. 
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The anomalous dimension due to SU(2) gauge boson 
exchange is obtained similalrly. If both fermions are dou- 



blets, the operator has the form C\ t A d 
and the anomalous dimension matrix is 



t A + C 2 1 ® 1, 
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where the (Ci,C 2 ) pairs are (Cqquj^, CQQ 2 iji), 

[CQQl2.fi, CQQ22,Ji), (CQLl.fi, CQL2,fi), 

(CLQl,fi,CLQ2,fi), and (CLL\,fi,ChL2,fi)- 

If one of the fermions is a weak doublet and the other 
is a weak singlet, the operator has the structure C 1 (g> 1, 
and the SU (2) anomalous dimension for C is identical to 
the Sudakov form-factor case, 



(106) 



1SU(2) 

If both fermions are weak singlets, then 

lsu { 2) = 0. (107) 
B exchange gives the diagonal contribution 



lu(i) 



OLl 
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(r l 2 + r/)(4U-6) + 8y l r / L t/i , 



(108) 



where Yi and Yf are the hypercharges of the initial and 
final representations. Note that Y(Ur) = 2/3, Y(Dr) = 
-1/3 and Y(E R ) = -1. 

The Higgs wavefunction graphs give the diagonal con- 
tribution 



iH(Q (t) ) 



16tt 2 2' 



(109) 
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to an operator for each field, and 



SU(2) gauge bosons. The matching matrix is 



lH{t R ) 



9t 
16tt 2 



(110) 



for each tu = field. This term breaks the flavor sym- 
metry in the anomalous dimension. The total anomalous 
dimension is the sum, 

7 = 7ff+7c/(i) +7sc/(2) +lsu(3), (HI) 
and is used to run the operators from /i = s to fi ~ Mz- 

C. Matching at the low scale to SCET 7 

At a low scale fi of order Mz (or m t ) one matches 
from SCETew with dynamical gluons and electroweak 
bosons onto SCET 7 with dynamical gluons and pho- 
tons, by integrating out the W and Z bosons. The elec- 
troweak symmetry is broken in SCET 7 so the operators 
in Eq. (|95|) must now be decomposed into separate SU (2) 
component fields. 

We start by considering the case where all particles 
have mass much smaller than m t , i.e. for all particles 
except the t-quark. This includes all the operators in 
Eq. ((95)) except those that contain and f/jj . The 
photon and gluon graphs are the same in SCETew an d 
SCET 7 and do not contribute to the matching condition. 
The W contribution depends on whether the particles 
invovled are SU(2) doublets or singlets. For the case of 
two doublets, consider the operators 

C QQ12J i[t a ] L <g> [t a ] L + C QQ22Ji [l] L ® [1] L (112) 

with i 7^ t, f t, which are two of the terms in Eq. (|M|) . 
For definiteness, let / = c and i — u. These operators 
match onto a linear combination of 



012 = 


[cL4n t iCLS\[uL2l' 1 U L i\ 




<5 22 = 


[cL4lt,C L3 }[d' L2 ^d' L1 } 




032 - 






042 - 






052 = 


[s' L4 lnCL3}[uL2j' 1 d' L1 } 




62 = 


[cL4l f ,s' L3 ][d' L2 j' i U L1 } 


(H3) 



where the flavor label represents the SCET 7 fields W^£. 

The matching from Eq. (fTP2|l in SCET EW onto 
Eq. {TJ5D in SCE T 7 is computed as in Sec. IVFH As 
shown in Sec. IV El the matching can be written as the 
sum of the Sudakov form-factor S'-matrix elements, even 
though we are considering W exchange, Z exchange 
and 7 exchange separately, and not summing over all the 
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022 
032 

C42 

052 
062 



R 



CqQ\2,cu 
0QQ22 ,cu 



(114) 



where Ca are the coefficients of Oi 2 ■ At tree-level R is 

1 



(115) 



At one-loop 
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47r sin 2 6 W 
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The Z exchange contribution is 
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47r sin 8w cos 2 #w 
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M 2 t 

rs = (g 2 Lu +g 2 Ls )F g (-s, M 2 z ,f) + Ag Lu g Ls \og-^- log - 

M 2 t 

= {gld + gL)Fg (-a, Mj ,n 2 )+Ag Ld g La log —f log - 
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(120) 



where 

9lc = ghu 
ghs = ghd 

are the couplings to the Z. The total one-loop matching 

ia p(i) , 0(1) 
is H w + H z . 

The remaining two operators in Cqq , 

Cqq = C QQUJi [T A t a ] L ®[T A t a } L 

+C QQ2lJt [T A ] L ®[T A ] L (121) 

match onto 

011 = [c L 4T A lfi C L3 }[uL2T A ^ULl} 

G n = [c L4 T A W3 }[d' L2 T A rd' L1 ] 
031 = [s' L4 T A lfl s' L3 }[u L2 T A ^u L1 ] 
041 = [sm^AsR^V^] 

051 = [^"V^^V^l] 

061 - [cHT\4 3 ][4rV«M]. (122) 

Since exchange leaves the color indices unaffected, 
the matching matrix is identical to Eq. (|115p . (|116[) . 
pi8p . and the matching relation is given by Eq. 1|114[) 
with the replace ment C l 2 — > 0q, Cqq i2iCU — > Cqqh jCU . 
The results Eq. (|115[) . (| 1 1 6[) hold for all cases where 
both fermions are doublets. If the final quark doublet 
is replaced by a lepton doublet, the coupling constants 
in Eq. (|119[) have the obvious replacement g^ c — > gLv, 
gLs -* gLe with 

1 



ffie = -- + sin 



(123) 



and similarly if the initial doublet is a lepton doublet, or 
both doublets are lepton doublets. 

The second case is where one fermion is a doublet 
and the other is a singlet. As an example, consider 
CQU2ji[l]L ® with f — c and i = u. This matches 
onto a linear combination of 

012 = [cLil^C L z][u R 2^U R i] 



02 



(124) 



The matching matrix is 



012 
022 



= RC, 



QU2,cu 



(125) 



where Ci 2 are the coefficients of Oi 2 ■ At tree- level R is 
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( 119 ) At one-loop 
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where F g is given in Eq. (1117|) . The singlet fermion Z 
couplings are 







gnu 


= -^sm 2 9 




1 „ 


gRd 


= g sin 


gRe 


= sin 2 ■ 
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Equations (|125[) . (|126j) . (|127| apply to all cases where 
one fermion is weak doublet, and the other is a weak sin- 
glet, with the obvious replacement of the Z charges for 
lepton doublets. Since electroweak exchange does not af- 
fect the color indices, the same matching matrix applies, 
for example, to the transition from CQuiji\T A \L®\T A \ R 
to 

0ii - [cL4T A ^c L3 }[u R2 T A ^u m } 

021 = [s' L4 T A lfl s' L3 }[u R2 T A ^u R1 ]. (129) 

The last case is if both fermions are weak singlets — 
take Cuu2,fi[i-]R <8> [1]r as an example with f — c and 
i = u. The operator matches to 



= [cRAltiC R3 ][u R2 ^u R1 ] . 
The one-loop matching condition is C = (1 + R 



R z )Cuu2 J- a 



(130) 
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Again, the same matching 



the matching between Cumji[T A 



coefficient holds for 
R ® [T A ] R and 

[cmT A ^ tl c R 'i][u R 2T A ^u R i], and the equations hold 
with an obvious substitution of Z charges if the quarks 
are replaced by leptons. 

Anomalous dimensions in SCET-y 

Finally, one computes the anomalous dimension of the 
operator in SCET 7 between /i ~ Mz and some low scale 
fig, at which point one takes operator matrix elements to 
compute the desired observables. The matrix elements 
of the initial state SCET fields in the proton are the 
usual parton distribution functions, and the final state 
fields are used to construct jet observables. The scale [iq 
is chosen to minimize logarithms in the matrix element 
computation. For LHC jet observables, it is of order the 
typical invariant mass of a single jet. 

The QCD anomalous dimensions given in Eq. (|101D . 
(|103[) . (|113p continue to hold for the case of two quarks, 
one quark and one lepton, or two leptons, respectively. 
The QED anomalous dimension is 



47T 



1 / 2 

2 Oh 



<7§ 



4 log -j- -6 



+4 (qxqi + q 2 q 3 ) log 4 (qiq 3 + q 2 q±) log - 

s s 

(132) 

where qi_4 are the charges of the four fields, and qi+qj, = 
12 + 94- The initial particle charges are q\ and —q 2 , and 
the final particle charges are —q 3 and 94. 



D. t-quark Production 



In processes involving the t-quark, m 2 /Myyz terms 
must be included in the loop graphs, as discussed in 
CGKM2. At a scale //, we transition to an effective theory 
in which the t-quark is represented by an HQET field, 
whereas the light quarks are still represented by SCET 
fields. Since mt is not much larger than My/ and Mz, 
it is convenient to make this transition at the same time 
that the W, Z and Higgs bosons are integrated out of 
the theory in the transition from SCETew to SCET 7 . 
This method was used in CGKM2 for the Sudakov form- 
factor of the t-quark, and allows one to include the com- 
plete m 2 /M^ vz dependence in the matching computa- 
tion. Here we apply the same procedure to the oper- 
ators relevant for t-quark production — the operators in 
Eqs. (|94l95p which contain either or t R fields. The 
initial state in proton-proton collisions contains only light 
quarks, so we will only look at operators with top-quarks 
in the final state and light quarks in the initial state. 

The matching at Q and the anomalous dimension be- 
low Q are mass independent, and identical to those for 
light quarks. The mt dependent terms give an additive 
correction to the low-scale matching matrices R^) z °f 
the previous section. There are also contributions R^ g 



to low-scale matching from the photon and gluon, be- 
cause of the transition to an HQET field for the t-quark. 
The graphs in SCETew use a SCET field for the t-quark, 
and those in SCET 7 use a HQET field for the t-quark, 
so there is a matching correction even for massless gauge 
bosons, as computed in CGKM2. 

The Higgs only contributes through wavefunction 
renormalization in SCET. The matching contribution 
from the Higgs is 



H(t L ) 
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-lS^ [ l Fh{Mhm2t)+ l Fh{Mhm2t) 

+-a(h t ,h t ) + -a(z t , z t ) + c(h t , h t ) + c(z t , 

+c(w t ,0) - b(h t ,ht) + b(z t , z t )] 

H(t L ) - ~ {MMwM) + a(w t ,0)) 



1 Vt 



2 16tt ; 
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\v 
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(!• 



where the functions are tabulated in Appendix [B] For 
each t£, t R or field, one adds H(ti), H(t R ) or 
H(b' L ) to the matching matrix. For example the oper- 
ator tL7 M tL6 / L 7 M 6^ gets the Higgs matching contribution 
2H(t L ) + 2H(b' L ). 

The matching for operators containing and a light 
quark doublet, Eq. (|112p with f = t and i = u,c, is given 
by Eq. (|116[) . with the additional additive contribution 
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from the Higgs, 
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from the W, 
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from the Z. 
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from the photon, where qt = 2/3 is the i-quark charge, 
and 



from the gluon where 
Wi 
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f F (w t ,0) - -a(w t ,0) - -c(w t , 0) 



^2 = f F (0,w t )--a(0,w t ) 
Ui = g 2 Lt f F (z t ,z t ) -\g\ t a(z u z t ) 
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(#lt + c(« t , zt) + 9Ltgmb(z t7 z t ) 
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9RtfF(z t ,z t ) 
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-g% t a(z t ,z t ) 



(git + 9m) c( z t,z t ) + gLt9mb(z t , z t ) 



(139) 

and the functions fp, a, b and c are tabulated in 
Appendix [Bj The matching matrix multiplied by 
(CQQu,t q ,CQQ2i,tq) gives the coefficients C k i, k = 
1,...,6, I = 1,2 of the operators in SCET 7 listed in 
Eq. (H3"|) . (fT2"2"j) with c ^ t and s ^ b' for the final state 
quarks, and the initial state flavors replaced by the two 
members of the light quark doublet q, (u,d') — > (c, s'), 
or (u,d') -> (tt,d'). Note that C 5 i, C 52 , C 6 i, C 62 are rel- 
evant for single-top production. 



The gluon matching SR g ± 1 is diagonal in color space, 



lg ^ it 

and does not mix the 1 ® 1 and T A ® T A operators. This 
follows from Eq. ([4*6")) and the additive nature of the mass 
corrections to the amplitudes. 

The matching for operators containing and a light 
quark singlet are given by Eq. (|127p with the additional 
terms 
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with Xi t 2 and f/i given in Eq. (|139|) . 

The matching for operators containg t R and a light 
quark doublet are given by Eq. (|127p with the additional 
terms 
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The matching for operators containing t R and singlet 
light quarks is given by Eq. (|13ip with the additional 
terms 
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Anomalous Dimension in SCET-y 
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The anomalous dimension in SCET 7 after integrat- 
ing out the electroweak gauge bosons and switching to 
HQET for the top quarks, is given by gluon and photon 
exchange. For t pair production, particles 3 and 4 are 
HQET i-quarks, and 7 can be obtained from Eq. (|44]) . us- 
ing the heavy-heavy anomalous dimension (71 of CGKM2) 
for exchange between (3,4), the heavy- light anomalous 
dimension (72 of CGKM2) for exchange between (3,4) 
and (1, 2) and the light-light anomalous dimension (73 of 
CGKM2) for exchange between (1,2). This gives Eq. (j6"l"T) 
for the QCD part of the anomalous dimension, with 
a —* a s , and group theory factors replaced by their 
SU(3) values, 
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The QED anomalous dimension is 



47T 



gf 41og-=--6 +<tf Ulog-5-4 



-8qiq t log - 
u 



(144) 



where qt — 2/3 and qi — 2/3,-1/3 is the charge of the 
light quark. 

For single-top production from the operators O51, Oqi, 
O52, Oe2i there is only one heavy quark in the final state, 
and 
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and the QED anomalous dimension 
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where q u = 2/3 and qd = —1/3 are the charges of the 
up- type and down- type quarks, respectively. 

This completes the computation of radiative correc- 
tions in the standard model. The formula? derived in 
this section will be used for the numerical computations 
in the next section. The only case we have not treated is 
when both initial and final particles are top quarks. This 
can be obtained from the case we have analyzed, with 
a heavy quark in the final state, by also adding heavy 
quark corrections terms for the initial quark. 





FIG. 10: Graphs contributing to gg — > qq in the full theory. 




FIG. 11: Operator contributing to gg — > qq in the EFT. 



E. 



99 -> qq, 91 -> 91 and gq -> gq 



The computations in this paper have been restricted to 
those involving external matter fields. In top-quark pro- 
duction and in jet production, processes involving exter- 
nal gluons are also important. Consider, for definiteness, 
the case gg — > qq. At the high-scale Q, the tree-graphs 
which contribute to gg — > qq are shown in Fig. [101 In 
the EFT, one generates a local operator which involves 
the fields q, q, and two collinear-gluon field strength ten- 
sors, shown graphically in Fig. [TTJ The QCD corrections 
involve studying operators with gauge field strength ten- 
sors, and will be discussed elsewhere. The QCD cor- 
rections are known from existing fixed-order computa- 
tions 53]. The new feature discussed in this article is 
the electroweak correction. If we restrict ourselves to the 
electroweak corrections alone, then we can compute these 
using the results in CGKM2. The gluon field strength ten- 
sor is an electroweak singlet, and so the ggqq operator 
in the EFT is equivalent to the electroweak singlet cur- 
rents qj^PL^Rq studied in CGKM2 and the running and 
matching corrections in the effective theory are identical, 
with the identification — Q 2 — » s. Thus the total radia- 
tive corrections are given by combining the known QCD 
corrections, with the electroweak corrections for the cur- 
rent given in CGKMl, CGKM2. The other important par- 
ton subprocesses which contributes to dijet production 
are gq — > gq, gq — » gq, and gg — > gg. For gluon-quark or 
gluon-antiquark scattering, the EFT operator is a ggqq 
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operator, and the electroweak corrections are the same as 
those for the Sudakov form-factor, with Q 2 — ► —t. For gg 
scattering, there are no electroweak corrections to the or- 
der we are working, since the gluons do not couple to the 
electroweak gauge bosons, and the radiative corrections 
can be computed using the known QCD corrections. 



F. Squark Production 

The techniques developed in the previous sections can 
be used to calculate the radiative corrections in a theory 
involving scalar particles in the final state such as SUSY. 
To perform a high precision computation requires speci- 
fying a particular supersymmetric theory, and computing 
the matching conditions and radiative corrections using 
the given SUSY particle spectrum. This is beyond the 
scope of the present work. 

To estimate the size of electroweak Sudakov correc- 
tions in squark production, we will compute the SU (2) 
corrections in the toy model of Sec. IV Dl assuming the 
squark is a doublet, and a — > the weak interaction 
coupling constant. This gives the expected size of elec- 
troweak Sudakov corrections in squark production. 



VIII. NUMERICS 

The formulae for the EFT computation of standard 
model scattering processes have been given in Sec. I VIII 
As discussed in Sec. IIIII the anomalous dimensions are 
integrated using the two-loop /3-functions, and we also 
include the known two-loop QCD anomalous dimen- 
sions [HI in addition to the one-loop results of Sec. IVIII 
The corrections have a very small dependence on the 
Higgs mass (much less than 1%). In the numerics, we use 
a Higgs mass of 200 GeV. The EFT coefficients should 
be run down to a scale fio of order a typical jet invari- 
ant mass. We have chosen to use fio — 30 GeV. The 
electroweak corrections are insenstive to this scale, be- 
cause the only electroweak correction below Mz is due 
to photon exchange. 

The matching corrections at the high scale Q are about 
2%, and dominated by the QCD contribution. The low- 
scale matching due to integrating out the W and Z 
is about 2%. Both matching corrections are not very 
strongly dependent on Q. The largest corrections are 
from the anomalous dimension running. These correc- 
tions grow rapidly with energy. The one-loop QCD cor- 
rections are very large, and reduce the rate by factors of 
3-30 in the range y/l between 1 and 5 TeV. The two- 
loop QCD cusp anomalous dimension reduces the rate 
by about 10% at vi = 5 TeV. This is smaller than the 
electroweak corrections, but not negligible. The two- 
loop non-cusp QCD anomalous dimension (the B term 
in Eq. (fT"5")0 increases the rate by about 2%. We have 
included the QCD two-loop cusp and non-cusp terms in 
the numerical results. The two loop cusp anomalous di- 



mension has been shown to be proportional to the one- 
loop result [34], and we use their K factor to determine 
the two-loop cusp anomalous dimension (the A term in 
Eq. ()15|) U The two- loop non-cusp anomalous dimension 
was determined in CGKM2 by comparing the EFT result 
with the two-loop results of Jantzen and Smirnov [13] • 
The two-loop cusp anomalous dimension also determines 
the two-loop contribution to the soft anomalous dimen- 
sion matrix 75. The non-cusp contribution vanishes, 
since 75 depends on differences of anomalous dimensions. 
The three-loop QCD cusp anomalous dimension contri- 
bution [54[ is less than 0.1%, and can be omitted. The 
one-loop electroweak anomalous dimension corrections 
are significant, ranging from 5% at 1 TeV to around 30% 
at 5 TeV. Higher order electroweak corrections, such as 
the two-loop electroweak cusp anomalous dimension are 
smaller than 0.1%. The numerical results are accurate 
at the one-percent level, so that the error in LHC cross- 
sections is dominated by other uncertainties, such as in 
the parton distribution functions. 

The EFT analysis neglected power corrections of the 
form M 2 /s, M 2 /t and M 2 /u. The dominant power cor- 
rections arise from one-loop QCD graphs, so we use the 
estimate (a s M 2 /-K) x l/(s,t, u) since the graphs have a 
color factor of (roughly) 4Cfch s / (47r). To keep the power 

corrections below 1% requires \J s, \t\ , |u| to be larger 

than about 200 GeV for light-quark processes, where the 
largest M is Mz, and larger than about 350 GeV for 
processes involving the top-quark. Note that we have in- 
cluded all power corrections that depend on ratios such 
as mz/m t or Mz/Mh, and not expanded in these ra- 
tios. There are tree-level power corrections due to gauge 
boson mass effects in the propagators, e.g the s-channel 
propagator s — M§ is approximated as s. These triv- 
ial effects cancel in our results, because we normalize all 
amplitudes to their tree- level values. 

The LHC cross-sections are given by using the coeffi- 
cients computed earlier to compute the parton scatter- 
ing cross-sections, and then convoluting them with par- 
ton distribution functions. For processes involving four- 
quark operators, the effective interaction at the low scale 
is a linear combination of two color structures, 
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Color-averaging over initial particles and color summing 
over final particles lead to a contribution to the cross- 
section which is proportional to an effective coefficient 
C, with 
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For qq 
cross-section is 
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66, the parton scattering 
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where Cll, etc. are the coefficients of the LL, etc. op- 
erators. The * denote partonic variables. The subscript 
s, t is a reminder that one uses the coefficients as com- 
puted in Sec. IVIll with annihilation channel kinematics. 
From this, one can compute hadronic cross-sections. For 
example, the dijet invariant mass distribution from the 
partonic subprocess uu — > dd is given by 



For qq 1 scattering processes not involving identical par- 
ticles the cross-section is 
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where M 2 is the dijet invariant mass, Ex is the transverse 
energy of the jet, r = M 2 /s, ^/s = 14 TeV is the LHC 
center of mass energy, and the parton luminosity function 
Lij is defined by 



~dT 



dx 



fP{x)ff\r/x) + ff\x)fl 2 \r/x) 
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(1 2) 

where f> : are the distribution functions for parton i 
in beams 1 and 2. 9 For the LHC, both are proton dis- 
tribution functions. The 1 + Sij is the symmetry factor 
for identical partons in the initial state. For the case 
q = q', e.g. uu — ► uu, Eq. (|149[) still holds, and the co- 
efficients Cll.rr get contributions from both the direct 
and crossed graphs. 

For identical particles, e.g. uu — > uu, the partonic 
cross-section has the schematic form 
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There is the direct channel as well as the crossed-channel 
with s t. For LL and RR, the crossed-channel ampli- 
tudes have the same fermion chiralities as the direct chan- 
nel, and are included as C, which includes both s <-> t 
and the crossing matrix. One has to add the amplitudes 
in the two channels before squaring. For LR and RL, the 
crossed-diagrams do not interfere because the initial and 
final chiralities do not match, and one adds the probabil- 
ities. 



and for qq' scattering 

da 1 
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The subscripts t, u, etc. are a reminder the one has to use 
the amplitudes of Sec. I VIII with the replacements s — > t, 
t — > u, etc. 

For identical quark scattering, qq — > qq, e.g. uu — > mm, 
the cross-section is 
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The 1/2 is from final state phase space for identical par- 
ticles. The initial state 1/2 is included in the parton 
luminosity function. 

There are 72 four-fermion amplitudes that have been 
computed in Sec. I VIII in the s-channel, not including 
those which are identical by flavor symmetry, and an- 
other 72 amplitudes in the i-channel, and we cannot 
plot them all here. We will choose some representative 
examples to illustrate the size of the radiative correc- 
tions in high energy LHC processes. Rather than plot 
the hadronic cross-sections, which involve convolutions 
over rapidly falling parton luminosities, we have chosen 
to plot the ratio of the partonic cross-sections da/dt to 
their tree-level values. From Eq. (|150p . it follows that 
this also gives the ratio of the hadronic cross-section to 
its tree-level value. We will also neglect the CKM ma- 
trix in the plots, since the flavor dependence of the elec- 
troweak corrections is small, and the CKM factors enter 
as off-diagonal CKM matrix elements multiplied by the 
difference of electroweak corrections between d! and s', 
and d! and b 1 . 



Plots 



2E T = V$ sin 8 and t = —s sin 8/2, where 8 is the center of 
mass scattering angle. Thus in Eq. {150}, a given Et values gets 
contributions from two values of t, or equivalently, one should 
symmetrize da/dt under t <-> u. We will plot dcr/df before sym- 
metrizing. 



Figure [T2l show the ratio der/di to its tree-level value 
for uu — * uu, cc, ti and bb as a function of for 

90° scattering, t = —s/2, including QCD and electroweak 
corrections. The radiative corrections are enormous, and 
reduce the cross-sections by 1.15-2 at = 500 GeV 
to a factor of 7-38 at ^/I — 5 TeV compared to the 
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FIG. 12: Rates for uu — » (dotted black), uu — > uu 

(solid cyan), uu — > cc (long-dashed red), uu — > ti (short- 
dashed blue), uu — » dd (dot-dashed green) and uu — > bb 
(double-dot-dashed magenta) as a function of y/§ in GeV at 
9 = 90°, normalized to their tree-level values without any 
electroweak corrections. Note the logarithmic scale. 



tree-level value, depending on the process. 10 The bulk of 
the correction is due to QCD effects. Some of the QCD 
corrections are included in parton shower Monte-Carlos, 
because gluon radiation from tree-level branching is re- 
lated to the LL Sudakov series. However, the electroweak 
corrections, and some of the QCD corrections are not 
included in the shower algorithms, so the Monte-Carlo 
results can have substantial (~ 50%) corrections. 

The uu rate differs from cc because of the crossed- 
channcl graph for identical particles. The difference be- 
tween cc and it, and between ss (not shown) and bb is 
due to top-quark mass effects. The uu — > [i + fi~ rate has 
smaller QCD corrections, since the final state is a color 
singlet. The anomalous dimension 7 is proportional to 
4L S — 6. At large values of s, the 4L S term dominates, and 
produces the large Sudakov (double-log) suppression. At 
smaller values of I, the —6 can compensate the 4L S term, 
leading to an enhancement of the cross-section. This 
leads to a flattening of the curves at the smallest val- 
ues of s. The cross-sections will decrease slightly if we 
continue the plot to even smaller values of s. This effect 
can also be seen in the plots of Ref . [9j . Figure Q2] show 
the radiative corrections to the angular distribution for 
uu — > uu, cc, ti and bb at vi = 1 TeV. There is 

about a factor of two variation in the radiative correction 
over the range —0.8 < t/s < —0.2. 

Figure [14] and Fig. [15] show the ratios da(uu — > 
ti)/dt/d&(uu — > cc)/dt and da(uu — > bb)/dt/d&(uu — » 
ss)/di as a function of y§ for different values of t (i.e. 



FIG. 13: Rates for uu — > (dotted black), uu — > uu 

(solid cyan), uu — > cc (long-dashed red), uu — ► ti (short- 
dashed blue), uu — > dd (dot-dashed green) and uu — > bb 
(double-dot-dashed magenta) as a function of i/s for Vs = 
1 TeV, normalized to their tree- level values without any elec- 
troweak corrections. 
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FIG. 14: The ratio (uu -> tt)/(uu -> cc) at t = -0.2s, (dotted 
blue), t — —0.35s (long-dashed red), t = —0.5s (solid black), 
t = —0.65s (double-dot-dashed magenta) and t = —0.8s (dot- 
dashed cyan) as a function of Vs in GeV. 



the scattering angle 9), including QCD and electroweak 
corrections. These ratios are unity in the absence of 
top-quark mass effects. There is a ~ 40% increase in 
the ti rate due to the top-quark mass. About —4% is 
from the Higgs contribution, —2% from mass effects in 
the low-scale electroweak matching, and the rest from 
mass effects in the QCD matching at m t and running 
below mt- There is a much smaller enhancement of 
da(uu — * bb)/dt/da(uu — * ss)/dt due to virtual top- 
quark effects in the ^-sector. 11 Figures. [TBI and [T7l show 



10 Note that the parton luminosity is falling by about four orders 
of magnitude over the same range. 



11 Even though our individual radiative corrections have corrections 
under 1%, ratios such as da(uu — > bb) /dt /da(uu —* ss)/dt have 
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FIG. 15: The ratio (mm -> bb)/(uu -> ss) at i = -0.2s, 
(dotted blue), i = —0.35s (long-dashed red), t — —0.5s 
(solid black), t — —0.65s (double-dot-dashed magenta) and 
i = —0.8s (dot-dashed cyan) as a function of \f§ in GeV. 



FIG. 17: The ratio of (dd -> bb)/(dd -> ss) at f = -0.2s, 
(dotted blue), t = —0.35s (long-dashed red), t — —0.5s 
(solid black), t = —0.65s (double-dot-dashed magenta) and 
i — —0.8s (dot-dashed cyan) as a function of \/I in GeV. 
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FIG. 16: The ratio (dd -> tt)/(dd -> cc) at 1= -0.2s, (dotted 
blue), t = —0.35s (long-dashed red), t = —0.5s (solid black), 
t = —0.65s (double-dot-dashed magenta) and t = —0.8s (dot- 
dashed cyan) as a function of yS in GeV. 



the corresponding results for dd — > ti, bb. 

The plots discussed above include QCD and elec- 
troweak corrections. To show the importance of elec- 
troweak corrections, we show in Fig. 1181 the same pro- 
cesses as in Fig. but instead of plotting the ratio of 
the partonic cross-section to the tree-level value, we plot 
the ratio of the cross-section to the value including only 
QCD corrections, i.e. with ai.2 — > 0. 12 This ratio shows 
the additional effect of electroweak corrections beyond 



much smaller errors, so that the deviation from unity in Fig. 1151 
is a real effect. 

In qq — > fi + fi~ , we include tree-level electroweak exchange, and 
keep a s ou,2 terms in the one-loop matching, but drop order af 2 
terms. 



1.1 
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FIG. 18: Electroweak corrections to mm — > /i + fi~ (dotted 
black), mm — » mm (solid cyan), mm — > cc (long-dashed red), 
mm — » ti (short-dashed blue), mm — > dd (dot-dashed green) 
and uu — ► bb (double-dot-dashed magenta) as a function of 
vl in GeV at 9 — 90°. The large corrections for mm — > dd 
arise from the t-channel W exchange graph. 



the QCD corrections, which have been computed pre- 
viously. The electroweak corrections are significant, in- 
creasing from (-4)-(-22)% at 1 TeV to (-18)-(-32)% 
at 5 TeV, depending on the process. The electroweak cor- 
rections to the angular distribution are shown in Fig. 1191 
There are 10-30% variations in the corrections in the 
range -0.8 < i/s < -0.2 for ^fs = 1 TeV. 

The electroweak corrections alone (defined as just dis- 
cussed) for lepton pair production from w-quark and d- 
quark annihilation are shown in Fig. [201 and Fig. [5T] 
for different values of i. At \fl = 1 TeV, the corrections 
range from (0.4)-(-14)%, increasing to (-13)-(-32)% 
at vl = 5 TeV. The electroweak corrections also change 
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FIG. 19: Electroweak corrections to uu — > (dotted 
black), uu — * uu (solid cyan), uu —> cc (long-dashed red), 
uu — > ti (short-dashed blue), uu — » dd (dot-dashed green) 
and uu — > bb (double-dot-dashed magenta) as a function of 
i/s for vi = 1 TeV. The large corrections for uu — > dd arise 
from the f-channel W exchange graph. 
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FIG. 20: Electroweak corrections to uu — ^ fi + fi~ at t = 
—0.2s, (dotted blue), t = —0.35s (long-dashed red), i — —0.5s 
(solid black), t = —0.65s (double-dot-dashed magenta) and 
t = —0.8s (dot-dashed cyan) as a function of y/1 in GeV. 



the angular distribution of the lepton pairs. Fig. [2U and 
Fig. [23] show the t dependence of the cross-section for 
different values of s. The angular dependence is approxi- 
mately independent of s. The reason is that the dominant 
t dependence arises from the soft anomalous dimension 
75, which is a function only of the dimensionless ratio 
i/s. The angular dependence of the electroweak correc- 
tions differ for uu — > an d dd — > . 

The electroweak corrections to heavy quark production 
via u and d quark annihilation are shown in Figs. [Ml 
Fig. [25l Fig. [26l Fig. [27] for i-quark production, and 
Fig. H3 Fig. [H Fig. (2D] Fig. ELD for 6-quark production. 

Electroweak corrections to heavy quark pro duction 
have been computed previously [55l [5rj l57l. l58l. l59l Itiol. [r3ll. 



FIG. 21: Electroweak corrections to dd — » /-i + fi~ at I = —0.2s, 
(dotted blue), t = —0.35s (long-dashed red), £ = —0.5s 
(solid black), t = —0.65s (double-dot-dashed magenta) and 
i — —0.8s (dot-dashed cyan) as a function of in GeV. 




FIG. 22: Electroweak corrections to uu — > at vi 

1 TeV, (dotted blue), V§ = 2.5 TeV (long-dashed red) 
andy/I = 5 TeV (solid black) as a function of i/s. 




FIG. 23: Electroweak corrections to dd — > fl ft at Vs = 
1 TeV, (dotted blue), v 7 ! = 2.5 TeV (long-dashed red) 
andv^ = 5 TeV (solid black) as a function of t/s. 
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FIG. 24: Electroweak corrections to uu — > ti at t = —0.2s, 
(dotted blue), i = —0.35s (long-dashed red), t = —0.5s 
(solid black), t = —0.65s (double-dot-dashed magenta) and 
i = —0.8s (dot-dashed cyan) as a function of \f§ in GeV. 
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FIG. 25: Electroweak corrections to dd — * ti at t = —0.2s, 
(dotted blue), t = —0.35s (long-dashed red), t = —0.5s 
(solid black), t — —0.65s (double-dot-dashed magenta) and 
i = —0.8s (dot-dashed cyan) as a function of y/1 in GeV. 
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FIG. 27: Electroweak corrections to dd — > ti at vs = 1 TeV, 
(dotted blue), yS = 2.5 TeV (long-dashed red) andvl = 
5 TeV (solid black) as a function of t/s. 
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FIG. 28: Electroweak corrections to uu — > bb at t = —0.2s, 
(dotted blue), t = —0.35s (long-dashed red), t — —0.5s 
(solid black), t = —0.65s (double-dot-dashed magenta) and 
t = —0.8s (dot-dashed cyan) as a function of \f§ in GeV. 
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FIG. 26: Electroweak corrections to uu — * ti at Vs = 1 TeV, 
(dotted blue), %/I = 2.5 TeV (long-dashed red) &iid\fi — 
5 TeV (solid black) as a function of t/s. 



FIG. 29: Electroweak corrections to dd — > bb at t — —0.2s, 
(dotted blue), t = —0.35s (long-dashed red), t — —0.5s 
(solid black), t — —0.65s (double-dot-dashed magenta) and 
t — —0.8s (dot-dashed cyan) as a function of \/l in GeV. 
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FIG. 30: Electroweak corrections to it it — > bb at v* = 1 TeV, 
(dotted blue), %/I = 2.5 TeV (long-dashed red) andvl = 
5 TeV (solid black) as a function of t/s. 
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FIG. 31: Electroweak corrections to dd — > bb at \/3 = 1 TeV, 
(dotted blue), %/I = 2.5 TeV (long-dashed red) and\/J = 
5 TeV (solid black) as a function of t/s. 



(32, 63]. We find the same qualitative behavior the elec- 
troweak corrections give a small (~ —6%) suppression, 
and the QCD corrections give a large (~ 50%) enhance- 
ment. 

The above plots have been for s-channel proceeses. 
There are also t-channel parton subprocesses that con- 
tribute to dijet production. Rather than go through these 
in detail, we show two illustrative plots: Fig.l32lshows the 
electroweak corrections to uu — > uu, ud — > ud, dd — > dd 
and ud — > ud (which is equal to du — > du) as a function of 
\/I for 90° scattering, and Fig. [33] shows the angular de- 
pendence of the electroweak corrections at = 1 TeV. 

There are also scattering processes involving external 
gluons. For gg — > gg, gq — ► gg, and gg — ► gg, we have 
only computed the electroweak part of the correction, 
which is equal to that for the Sudakov form-factor. For 
the s-channel procees gg — > qq, the electroweak correc- 
tion only depends on \/§, and is shown in Fig. 1341 The 
same plot also gives the electroweak correction to the t- 
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FIG. 32: Electroweak corrections to uu — > nit (dotted black), 
ud — > ud (long-dashed red), dd — > dd (solid blue) and ud — > 
ud, du —> du (dot-dashed green) as a function of y/1 in GeV 
at 9 = 90° . 
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FIG. 33: Electroweak corrections to uu — > uu (dotted black), 
ud — > ud (long-dashed red), dd — > titi (solid blue) and ud — > 
wd, dit — > du (dot-dashed green) as a function t/s at = 
1 TeV. 



channel scattering processes gq — > and gg — > gq as a 
function of \f—i, by crossing symmetry. The imaginary 
parts from the logarithmic branch cuts in the s-channel 
amplitude do not change the absolute value of the am- 
plitude. 

Finally, we show the electroweak corrections for squark 
production. As discussed earlier, we use the electroweak 
correction for squark production in the toy theory, with 
the gauge coupling constant set equal to ai of the stan- 
dard model. This gives an indication of the size of elec- 
troweak corrections to squark production in supersym- 
metric extensions of the standard model. A more pre- 
cise computation depends on the specific scenario. The 
electroweak correction to squark production is shown in 
Fig. [35] for a squark mass of 250 GeV. The radiative cor- 
rection to the angular distribution is much smaller than 
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FIG. 34: Electroweak corrections to gg — > uu, cc (dotted 
black), gg — > dd,ss (long-dashed red), gg — > ti (solid blue) 
and gg — > bb (dot-dashed green) as a function of Vs in GeV. 
The electroweak corrections are independent of I to the or- 
der we are working. The same plot also gives the electroweak 
corrections to gq —* gq and gq — > gq, as a function of sf—i. 
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FIG. 35: Electroweak corrections to qq — > qq* in the toy 
theory at t — —0.2s, (dotted blue), t = —0.35s (long-dashed 
red), i = —0.5s (solid black) as a function of ■s/S in GeV. The 
rate is symmetric under 8 — > 180° — 6. 



for fermions. For discovering squarks, the only correction 
which matters is that at threshold, = 2m g - since the 
parton luminosity falls steeply with s. The electroweak 
corrections give a small (5%) enhancement of the rate for 
s near threshold. 
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FIG. 36: Electroweak corrections to qq — > qq' in the toy 
theory at Vs = 1 TeV as a function of —t/s at y§ — 1 TeV 
(black), sf§ = 2.5 TeV (long-dashed red) and Vs = 5 TeV 
(dotted blue). 



the Sudakov logarithms have been summed using renor- 
malization group methods. The EFT also properly sums 
mixed higher order logarithms that depend on both a s 
and Oi\ t 2, as well as those that depend on the top-quark 
Yukawa coupling. We have checked that our results agree 
with previous results when expanded in powers of a. 

The electroweak corrections can be important for LHC 
processes, particularly in searches for new physics that 
look for deviations from the standard model. The cor- 
rections vary in size from about (0.4)-(— 14)% at 1 TeV to 
about (— 13)-(— 32)% at 5 TeV, and need to be included 
to obtain LHC cross-sections with accuracies under 10%. 

We have also shown that the radiative corrections to 
four-quark operators are given in terms of those for two- 
quark operators by summing over pairs of particles. The 
relation between this and factorization, and with the two 
loop soft anomalous dimension of Aybat et al. [34[ was 
discussed in Sec. I VII Further work on this important 
topic is in progress. 
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APPENDIX A: MATCHING AT Q INCLUDING 
THE CASE OF IDENTICAL PARTICLES 



IX. CONCLUSIONS 

This paper extends the analysis of two previous pub- 
lications P, , and gives detailed numerical results for 
radiative corrections to high energy scattering processes 
in the standard model. The electroweak and QCD cor- 
rections have been computed using EFT methods, and 



In this appendix, we summarize the matching compu- 
tation at scale Q, including the case of identical particles. 
We start with the case of an SU(N) gauge theory with 
left-handed fermions in the fundamental representation. 

The tree-graph in Fig. [37] gives 

A« = — [u 4 ^T a v 3 ) L [v 2l »T a Ul ] L (Al) 
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FIG. 37: 



which is written as 

A (s) = 4™ 
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A LL 



All (T a (g>T a ) 



(A2) 



factoring out the color structure from the spinor structure 
of the graph. The left T a is contracted with the color 
indices of particles 4 and 3, and the right T a with those 
of particles 2 and 1 . 

If the initial and final particles are identical, there is 
also the t-channel graph in Fig. 1381 which gives 
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FIG. 38: 

[u^T a Ul ] L [v 2 j»T a v 3 ] L 
[u A ^u l ] L [v 2 ^v 3 ] L (T a ®T a ) c (A3) 



and the relative minus sign is from Wick's theorem. The 
subscript c indicates that the color structure is in the 
crossed i-channel. The left T a is contracted with the 
color indices of particles 4 and 1, and the right T a with 
those of particles 2 and 3. 

It is convenient to convert the i-channel graph to the 
standard basis used in the paper. The i-channel color 
structure can be converted to the s-channel using the 
SU(N) color crossing matrix 

[ {T a <g> T a ) c (1 ® l) c ] = [ (T a <g T a ) (1 ® 1) ] M N . 

(A4) 



The color Fierz identity 
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can be written as 

(T a <g> T a ) 
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Using this and the same equation with direct and crossed 
channels exchanged, one finds 
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with Af^r = 1. There is no color crossing matrix required 
for a 17(1) gauge theory. For 57/(2) and SU(3), 
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The spinor Fierz is 



\U A ^U X \ L [V27 AI «3] L 
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so Eq. 
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Comparing with Eq. (|25|) in Sec. IVJ we see that the 
s-channel contribution to the matching coefficient is 







and the t-channel contribution is 
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The total contribution is 
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where the i-channel pieces should only be included for 
identical particles. 

This sets up the notation and procedure to be used for 
the one-loop matching computation. The full theory dia- 
grams of Fig. [5] were computed in order to match the full 
gauge theory onto SCET at n — Q. Dimensional regular- 
ization was used to regulate both the infrared and ultra- 
violet divergences, which are distinguished by subscripts 
on 1/e. The diagrams are computed with all masses set 
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to zero. The logarithms are expressed using the short 
hand notation 

L x = log — T , L x/y = log - (A14) 
ft* y 

for x,y = s, t, u. 

The first two vertex graphs of Fig. [6] each give a con- 
tribution of 



a 



All ( Cf — t.Ca 



1 



-— L s -L 2 + 3L s -8 + ^- 
em 6 



£ UV e iR e IR 
2 " 



(A15) 



The next two vertex graphs in Fig. [6] each involve a triple 
gauge boson coupling, and give 



V a 



a A ° A 
- AlL ~ 



2 + L, 



euv eiR 



(A16) 



The s-channel box graph in Fig. (|39p with all fermions 
left-handed gives 



FIG. 39: 
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V b = a 2 1 2 {s, t) All 
= a 2 I 2 (s,t)A LL 

where 



MM) = MM)--/(M) 
s 
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iV 2 - 1 iV 2 -2 
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(A17) 
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FIG. 40: 



The gauge boson self-energy graphs combine to give a 
contribution of 



V s = —All\C a 



31 5. 



3euv 9 3 



+T F riF 
+T F n s 



4 20 4, 
h -U 

3euv 9 3 

1 8 1, 



3euv 9 3 
and the wavefunction graph is 

« 2 , r 1 1 

V w = —C f All -= + = — 

s 2euv 2eiR 



(A20) 



(A21) 



The sum of all of the diagrams of Fig [5] including the 
gauge boson self-energy graphs and the wavefunction 
graphs is 

VWi = 2V v + 2Vg + V b + V c + V s +W w (A22) 
and gives 

>,2 f 



LL 



"i 

2C F 



£ir £ir £ir 
2 



-U-L? + 3L,- 
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20 4, 

¥ + 3 l 



1. 



s(s + 2t) 
(s + t) 2 



L t/s 



2s 
s+1 



-t/s 



(A18) 



— L t/u -4L s L t/u -/(s,t) 
em 



The s-channel crossed-box graph in Fig. (|40p with all 
fermions left-handed gives 



.4 



2LL 



a 
s 



— L t/tt -4L s L t/u -/(s,t) 
eiR 
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( C d + C A ) 
4 

Ci 



(A23) 



V c = -a 2 h{s, u) All 
= -a 2 h{s,u)A L L 



C x 1® l + ~(C d -C A )T a ®T a 



N — 1 1 



which are the coefficients of A L L(T a <g> T a ) and ® 
1), respectively. 

The counterterms of the full theory have been used to 
cancel to the 1/euv terms and the remaining poles are 
(A19) all 1/eiR infrared divergences. These infrared divergent 
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terms agree with the ultraviolet divergences in the effec- 
tive theory. The finite parts of Eq. (|A23| give the high 
scale matching condition at p ~ Q. 

Equation (|X23|) gives the one-loop matching result for 
fcrmions which are distinguishable. If the fcrmions are 
identical, then there are also one-loop graphs in the 
crossed channel, analgous to the crossed channel tree 
graph Fig. [351 They are obtained by the crossing relation 
Eq. (|A12[) used earlier for the tree- level graphs. 

The one-loop matching conditions for initial and fi- 
nal fcrmions of the same chirality (i.e. LL or RR) is 
Eq. (j A23|) . If the fcrmions have opposite chirality, then 
one can obtain the matching coefficients using charge 
conjugation. The right handed field ipR is replaced by 
the charge conjugate field ^>£. This reverses the sign of 
the fermion arrow on the fermion line, and exchanges the 
box and crossed box graphs. One can now use Eq. (|A23[) 
for the same-chirality case, and then use charge conju- 
gation on the final amplitude to rewrite the ipl spinors 
in terms of the original ipn spinors. The result of this 
procedure is that the matching Eq. (|A23[) for opposite 
chirality (i.e. LL or RR) is given by Eq. (|A23[) with the 
replacement C\ — » —C\, Ca — ► — Cd and t <-> u. 



APPENDIX B: PARAMETER INTEGRALS 

The parameter integrals tablulated below arise from 
vertex and wavefunction graphs where the gauge boson 
has mass M, the external particle has mass m ex t, and 
the internal particle has mass m. They depend on the 
variables w — m 2 xt /M 2 and z = m 2 /M 2 . For any func- 
tion f(w,z) defined below, we define the corresponding 
function of a single argument by f(z) = f{z,z). In 
the standard model, where the only fermion with mass 
comparable to the gauge boson masses is the top quark, 
we need the integrals f(z,z), f(z,0) and f(0,z), with 
z = m 2 /M 2 v ,m 2 /M 2 ,m 2 /M% . 

For 4z > 1, the f(z,z) results can be ana- 
lytically continued using y/l — 4z — > iy/4z — 1 and 
tanh _1 ("v/l — 4z) — ► itan~ 1 (\/4z — 1). In each integral, 
the factors of i cancel, and the function remains real. 



The f(w, 0) formulae are given by using f(w + i0 + , 0) for 
w > 1. They have an imaginary part for large values of 



Fermions 



The gauge boson vertex graph leads to the integral 



f F (w,z) = 2 



, 1 — x , / 1 — x + zx — wx(l — x) 
ax log 1 



f F (z,z) = 2 



1 - x 



| - 2 J log(z) + ilog 2 (z) 



(Bl) 



+ 



— tanh -1 Vl - 4z 



-2 (tanh -1 y/l - Az) . (B2) 
The function /f(-z) = /f{z, z) was used in CGKM2. 



f F (0,z) = — 



TT 2 2Z 



3 1-z 



log*-2Li 2 (l-z) (B3) 



1 — w 

f F (w,0) = 2 + 2 log(l-w)-2Li 2 (w)(B4) 

w 

The inverse propagator including the gauge boson 
wavefunction graph is 

S- 1 = f [1 + A(p 2 )] - m oxt [1 + B(p 2 )] . (B5) 

The parameter integrals required are 

a{p 2 /M 2 ,m 2 /M 2 ) = A(p 2 /M 2 , m 2 /M 2 ) 



b(p 2 /M 2 ,m 2 /M 2 ) = P 2 ^(p 2 /M 2 ,m 2 /M 2 ) 

op z 



c{p 2 /M 2 ,m 2 /M 2 ) = p 2 ^( p 2 /M 2 ,m 2 /M 2 ) (B6) 

op z 

where m is the mass of the internal fermion, and the 
integrals are evaluated on-shell, with p 2 — m 2 xt , where 
TO 0X t is the mass of the external fermion. 

f 1 fl-x + zx-wx(l-x) 

a(w,z) = —2 1 da;(l — x) log 

Jo V 1 - x 

1 AJwz x(l — x) 
ax- t - 

,, 1 — x + zx — wx(l — X) 

1 , 2wx(l - x) 2 
ax -. 

(1 1 — x + zx — wx(l — X) 



b(w, z) 
c{w, z) 



(B7) 



a (z,z) = 
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i v / T^ 7 4i 
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1 - 4z + 2z 2 
2^2 



log z 



a(0, z) 
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logz 
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6(0, z) 
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log(l - w) 



, 4(32-1) , , . 

-4 + \ ' tanh -1 Vl-4z 
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, \ 2 o 2 1-5z + 5z 2 , i , — 

c(z,z) = --3+-: ; tanh" 1 Vl - 4z 

z z 2 Vl - 4z 

1 -3z + z 2 , 
H 5 log z 



c(0,z) 
c(w,0) 

The function 





2 , 2(1-11;), „ 
- - 1 + V 9 ; log 1 - to) 
to 



(BIO) 



/i F (z) = a(z,z)-2fc(z,z) + 2c(z,z) (Bll) 

was used in CGKM2 and is the wavefunction correction in 
a vector-like theory. 

The corresponding functions for radiative corrections 
due to a virtual scalar are 



a(w,z) = -a(w,z) 



b(w,z) = --b(w,z) 



c(w,z) = -c(w,z) 



(B12) 



and 



h F {z) = a(z,z) -2b(z,z) + 2c(z,z) (B13) 
was used in CGKM2. 

2. Scalars 

The gauge boson vertex graph for scalar particles leads 
to the integral 



fs(w,z) 
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Jo 



dx 



(2 — x) 1 — x + zx — wx(l — x) 



■log 
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= I - ( 1 - log(z) + \ log 2 z 



fs(0,z) = 



+ V1 4z tanh- 1 ( A /T~4i) 

z 

-2 (tanh -1 Vl -Azf 

2 

n L- + -J—\ogz-2U 2 (l-z) 

6 1 — Z 



1 — w 

fs{w,0) = 1 + — — log(l-w)-2Li 2 H 
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Scalar wavefunction renormalization due to gauge bo- 
son exchange gives the integral 



Iisitr. j n.rl (3x 2 - 6x + 4) 
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Scalar wavefunction renormalization due to scalar ex- 
change gives: 



h s (w,z) 



dx 
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APPENDIX C: ERRATUM 

The low-scale matching for the i-quark in CGKM2 is 
incorrect. The corrected expressions are 



[^l 2 w n ]^p R [wUf tP1 ] 



ait V2 ^ P R t Vl , 



(CI) 
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where the matching coefficients ai_3 are given by 

logaiK) = - . 2 " cm [gl t F g (Q, M z , m t ) + 2U,] + °^ (±) [F g (Q, M Wt m t ) + 2W{\ 

47r sin Ow cos Ow Air sm B w \ l J 

( a s 4 a cm 4\ / n 2 \ „ , , 

loga 2 (m t ) = - . " cm gl b F g (Q,M z ,m t ) + - ^ [F ff (Q, M Wl m t ) + 2W 2 ] + 2ff(6' L ), 

47r sm t^vF cos 2 #w 47r sm \ 2 / 

and the required functions Ui 2, A"i 2, HUl), H(tn), and H{b' L ) are given in Eqs. (|117p . (|133p . (|139p . 
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